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ABSTRACT 

Aims. We describe various expansion schemes that can be used to study gravitational clustering. Obtained from the 
equations of motion or their path-integral formulation, they provide several perturbative expansions that are organized 
in a different fashion or involve different partial resummations. We focus on the two-point and three-point correlation 
functions, but these methods also apply to all higher-order correlation and response functions. 

Methods. We present the general formalism, which holds for the gravitational dynamics and for any similar models, 
such as the Zeldovich dynamics, that obey similar hydrodynamical equations of motion with a quadratic nonlinearity. 
We give our explicit analytical results up to one-loop order for the simpler Zeldovich dynamics. For the gravitational 
dynamics, we compare our one-loop numerical results with numerical simulations. 

Results. We check that the standard perturbation theory is recovered from the path integral by expanding over Feynman's 
diagrams. However, the latter expansion is organized in a different fashion, and it contains some UV divergences that 
cancel out as we sum all diagrams of a given order. Resummation schemes modify the scaling of tree and one-loop 
diagrams, which exhibit the same scaling over the linear power spectrum (contrary to the standard expansion). However, 
they do not significantly improve over standard perturbation theory for the bispectrum, unless one uses accurate two- 
point functions (e.g. a fit to the nonlinear power spectrum from simulations) . Extending the range of validity to smaller 
scales, to reach the range described by phenomenological models, seems to require at least two-loop diagrams. 

Key words, gravitation; cosmology: theory - large-scale structure of Universe 



1. Introduction 

According to standard cosmological scenarios, the large- 
scale structures of the present Universe arise from the am- 
plification by gravitational instability of small primordial 
fluctuations (Peebles 1980). Then, from observations of the 
recent Universe at rcdshifts z < 5, through galaxy sur- 
veys (Tegmark et al. 2006; Cole et al. 2005), weak-lensing 
studies (Massey et al. 2007; Munshi et al. 2007), or stud- 
ies of baryonic oscillations (Eiscnstcin et al. 1998, 2005), 
we can derive constraints on the cosmological parameters 
and on the properties of those primordial fluctuations. This 
requires accurate theoretical predictions for the evolution 
of gravitational clustering to compare cosmological models 
with observations. On large scales, where density fluctua- 
tions are small, linear theory is sufficient. On small scales, 
deep into the nonlinear regime associated with virialized 
structures such as galaxy clusters, one can use numerical 
simulations or phenomenological models, such as the halo 
model (Cooray & Sheth 2002), which allow for an accuracy 
near 10%. However, many observational probes (such as 
weak lensing or baryonic oscillations) are mostly sensitive 
to intermediate scales (wavenumbers k ~ Q.lh Mpc~^ at 
z ^ 1) that might be studied through analytical methods, 
such as perturbative expansions. 

Thus, there has recently been a renewed interest 
in building analytical methods that could improve over 
the standard perturbation theory and provide reliable 
predictions in this weakly nonlinear regime. Crocce & 
Scoccimarro (2006a, b) present a partial resummation of 



the diagrammatic series associated with the response func- 
tion, which recovers the expected decay into the nonlin- 
ear regime, in agreement with numerical simulations. This 
allows one to obtain a more accurate prediction for the 
matter power spectrum than with the usual perturbation 
theory (Crocce & Scoccimarro 2007). On the other hand, 
Valageas (2004, 2007a) described the path-integral formal- 
ism associated with the cosmological gravitational dynam- 
ics. This allows one to apply the usual tools of field theory, 
such as large-iV and semi-classical expansions. This pro- 
vides several methods to perform partial resummations of 
the usual perturbation theory. Next, Matarrese & Pietroni 
(2007a, b) proposed to use the dependence of the system 
on a high-/c cutoff to derive some resummation schemes 
through the path-integral formalism. An alternative ap- 
proach developed by McDonald (2007), based on a renor- 
malization group technique, also permits going beyond the 
standard perturbative expansion. 

These works have focused on two-point functions, such 
as the matter power spectrum or two-point correlation, 
since the latter are the main quantities obtained from ob- 
servations. However, the additional information included in 
higher-order statistics, such as the three-point correlation, 
can prove very useful in breaking degeneracies (Sefusatti 
et al. 2006; Munshi et al. 2007). Therefore, we describe 
in this paper how such methods can be used to compute 
higher-order functions. We focus on the bispectrum, which 
is the quantity of most practical interest, since higher-order 
statistics obtained from observational surveys are increas- 
ingly noisy. We consider both the exact gravitational dy- 
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namics and the simpler Zeldovich dynamics. Indeed, the 
equations of motion associated with both systems are simi- 
lar and these expansion schemes apply in the same manner 
to both (Valageas 2007b). Then, the simpler Zeldovich dy- 
namics can be used as a convenient benchmark to track 
possible errors in the derivation of the equations associated 
with various schemes. 

This article is organized as follows. First, in Sect. [2] we 
recall the equations of motion associated with cosmologi- 
cal gravitational clustering. We focus on weakly nonlinear 
scales hence we restrict ourselves to the hydrodynamical 
description. We write the path-integral formulation of this 
system and we briefly present two large- methods which 
can be derived within this formalism. We also recall the 
main properties of the simpler Zeldovich dynamics. Next, 
we present in a general manner various perturbative ex- 
pansions that can be applied to both cosmological dynam- 
ics (or any equation of motion of the same form, i.e. with 
a quadratic nonlinearity). Thus, we first recall the stan- 
dard perturbative expansion over powers of the linear den- 
sity field in Sect. [3l we describe the equivalent expansion 
that can be derived from the path-integral formalism by ex- 
panding over the interaction term of the action in Sect. SI 
we present the resummations associated with the large- 
methods in Sect.O and we briefly discuss in Sect. [S] the re- 
summation associated with the random advection of small- 
scale structures by the larger-scale coherent flow (involv- 
ing a high-fc limit). Then, we apply these methods to the 
Zeldovich dynamics in Sect. [71 up to one-loop order. Indeed, 
the expressions obtained in this case are simpler than for 
the gravitational dynamics. This allows us to show explic- 
itly the connection between the different methods and to 
point out the possible divergences which may appear in in- 
termediate steps as diagrams are gathered in different man- 
ners. Finally, we apply these methods to the actual gravita- 
tional dynamics in Sect.[8l again up to one- loop order, and 
we compare their predictions with numerical simulations. 
We conclude in Sect. [51 

The reader who is not directly interested in the theo- 
retical aspects of gravitational clustering (e.g. the discus- 
sion of the various formalisms and perturbative methods 
that may be built, and the problems, such as UV diver- 
gences, that one may encounter) might first read Sect. 18.41 
and Figs. 1^011^ There, we compare the standard pertur- 
bation theory with the most efficient resummation scheme 
considered in this paper (at one- loop order), as well as a 
recent phenomenological model and numerical simulations. 

2. Equations of motion 

2.1. Fluid dynamics 

In this section, we recall the equations that describe the dy- 
namics of gravitational clustering in a cosmological frame- 
work. As in Valageas (2007a), we focus on weakly nonlinear 
scales where we can use the single-stream approximation. 
Then, both cold dark matter and baryonic matter are de- 
scribed as a pressureless fluid governed by the continuity 
and Euler equations, supplemented by the Poisson equa- 
tion for the gravitational potential (f) (Peebles 1980), 



(3) 



|^ + V.[(l + 5)v]=0, 

T7^ + -t- (V.V)V = - 

OT 



where r ^ Jdt/a is the conformal time (and a the scale 
factor), Ti. = dlna/dr the conformal expansion rate, and 
f2in the matter density cosmological parameter. Here, S is 
the matter density contrast and v the peculiar velocity. 
Since the vorticity field decays within linear theory (Peebles 
1980), we take the velocity to be a potential field so that v 
is fully specified by its divergence 6 = V.v. It is convenient 
to work in Fourier space, which we normalize as 

dx 



^<5(x). 



(4) 



and to introduce the time coordinate rj defined from the 
linear growing mode _D+ as 

T]^\nD+{T) with D+{z^ 0)^1. (5) 

Then, as in Crocce & Scoccimarro (2006a), we define the 
two-component vector as 



■0i(k,?7) 
■02 (k, 7?) 

where /(r) is defined as 

dlnZ?! 1 dlnD, 



(6) 



(7) 



■^''^^ din a H dr 

Then, the equations of motion (P)-® read (Crocce & 
Scoccimarro 2006a; Valageas 2007a): 

0{x,x').ip{x') ^ Ksix;xi,X2).i>ixi)ip{x2), (8) 

where we have introduced the coordinate x — (k, 77, i) where 
i = 1, 2 is the discrete index of the two-component vectors. 
In Eq.® and in the following, we use the convention that 
repeated coordinates are integrated over as 

2 

Oix, x').^(x') = / dk'dry' J20u' (k, ry; k', 77')^,' (k', ??')• (9) 

i' — l 

The matrix O reads 



0{x,x')-. 




whereas the symmetric vertex Ks{x; xi,X2) — Ks{x; X2,xi) 
writes as 



Ks{x;xi,X2) 
with 

7i;i,2(ki,k2) 
7i';2.i(ki,k2) 

72:2.2 (kl,k2) 



Soiki -I- k2 - k)6D{rii 

X7f;ii,»2(kl,k2) 



■11)5 Dim - v) 



(ki 



k2).k2 



2fc| 
(ki -fk2).ki 

|ki+k2p(ki.k2) 



(11) 

(12) 
(13) 
(14) 



and zero otherwise (Crocce & Scoccimarro 2006a; Valageas 
2007a) . In the following, we use the common approximation 

- 1, (15) 



(1) 7 

which simplifies the analysis as the equation of motion ([8]) 

(2) no longer shows any explicit dependence on time rj. 
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2.2. Linear regime 

On large scales or at early times where the density and 
velocity fluctuations are small, one can linearize the equa- 
tion of motion ([5]) which yields 0.i(jl = 0. Then, the linear 
growing mode is 

M^)^e''SM(^iy (16) 

where SloO^) the linear density contrast today at redshift 
z = 0. As usual, we define the initial conditions by the linear 
growing mode (|16p and we assume Gaussian homogeneous 
and isotropic initial conditions defined by the linear power 
spectrum Plo(^): 

((5Lo(ki)<5Lo(k2)) = PLo{ki)SDiki + ks). (17) 

Then, the linear two-point correlation function CLixi, X2) 
reads as 

Cl{xi,X2) = {iPl{xi)'iPl{x2)) 

= e^'+^'PLo{ki)SD{ki + k2) J \ y (18) 

In previous papers (Valageas 2007a, b) we noted the two- 
point correlation as "G" , but since a restriction of the re- 
sponse function introduced below is also labeled "G" in the 
literature, we shall note correlation functions as "C" in this 
article. As in Valageas (2004), it is convenient to introduce 
the response function R{xi,X2)- Using the approximation 
(fTS]) it reads in the linear regime (Crocce & Scoccimarro 
2006a; Valageas 2007a) 

Rl{xi,X2) = 60(^.1 - k2) 9{r]i - 772) 

n^l32j+ 1 [-3 3 j|' (19) 

where the Heaviside factor 9{rii —772) ensures that causality 
is satisfied. Thanks to statistical homogeneity and isotropy 
the correlation functions G (both the nonlinear G as well 
as Gl) are symmetric and of the form 

C{xi,X2) = (5_D(ki +k2)Gii,i2(fci;77i,772), (20) 

with 

Ct^,t2{k;m,V2) = C^^,^,{k■,r|2,m), (21) 

whereas the response functions R (again both the nonlinear 
R as well as i?^) are of the form 

R{xi,X2) = (5_D(ki - k2)6'(r/i - ?72)i?Ji,j2 (^1; ^i. 'ya)- (22) 

Moreover, the linear two-point functions obey (Valageas 
2007a,b): 

0(x,z).GL(z,y) = 0, 0{x,z).RL{z,y) = dD{x-y). (23) 
This can also be checked from expressions ([T^. 



where we have introduced an imaginary auxiliary field A(x), 
and the action S[>p,X\ reads a43 

S[il;, A] = A.(O.V - K^.iji;) - ^A.A/.A, (25) 

where A/ is the two-point correlation of the initial condi- 
tions taken at time rji. This matrix disappears in the final 
equations when we take the limit 77/ —^ —00. Moreover, the 
first and second moments of the field A obey 

(A>-0, (AA)=0, i^Pixi) Xix2)) ^ RixuX2). (26) 

The action ([^5]) can serve as a basis for several approaches 
borrowed from field theory. For instance, it was used in 
Matarrese & Pietroni (2007a, b) to derive a resummation 
scheme based on the dependence of the system on a high- 
k cutoff (see also Valageas 2007b for a discussion of this 
approach in the case of the Zeldovich dynamics). Another 
approach relies on standard "large- A'^" expansion schemes 
(Valageas 2007a, b). Thus, generalizing the problem to A^ 
fields or simply multiplying the exponent in expression (j24p 
by a factor N, one can derive various expansions over pow- 
ers of 1/A^ (which agree up to the order 1/A^^ of the ex- 
pansion but differ by higher-order terms as they involve 
different partial resummations). This path-integral formal- 
ism can also be written for the integral form of the equation 
of motion, as discussed in Valageas (2007b). 

2.4. Direct steepest-descent method 

A first large- A^ expansion is provided by a direct steepest- 
descent method, which expands the path integral p4)l 
around a saddle-point as described in Valageas (2004, 
2007a). It introduces auxiliary correlation and response 
functions Go and Rq defined by evolution equations that 
happen to be identical to the linearized equations (j23p . 
Thus, we have Go — Cl and Rq — Rl- (This is specific 
to the hydrodynamic description, where {'ip) — 0, and does 
not hold if we use the Vlasov equation of motion to describe 
the system.) Then, the nonlinear correlation and response 
functions G and R can be obtained from the coupled linear 
equations 

0{x,z).C{z,y) = J:ix, z).C{z,y) + nix, z).R^{z,y) (27) 
0{x, z).R{z, y) = 6D{x-y) + Y.{x, z).R{z, y) (28) 

where we have introduced the self-energy matrices S and 11 
that can be expressed in terms of Go and Rq (hence in terms 
of Gl and Rl) - Equations (P7)l - (P51) are actually exact (they 
may be seen as the definition of the self-energy matrices) 
and the order of the large- A^ expansion only enters 

the expression of the self-energy, which can be obtained 
from a series of diagrams. At one-loop order, there is only 
one contribution which reads as 



2.3. Path-integral formalism and large-N expansions 

As described in Valageas (2007a, b), the statistical proper- 
ties of the density and velocity fields can be obtained from 
the generating functional 

Z[j] = (e^-^) ^ f [di^][dX] e^-^-si^^]^ (24) 



T.{x, y)^AKs{x] xi,X2)Ks{z; y, Z2)Rl{xi, z)Cl{x2,Z2) (29) 
n(a;, y)^2Ks{x; xi,x2)Ks{y; yi, y2)CL{xi, yi)CL{x2, y2) (30) 



Note that the action S only takes such a simple form because 
the Jacobian, associated with the change of variables from the 
initial condition tpL to the nonlinear field -i/;, can be shown to be 
an irrelevant constant. 
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2.5. The 2PI effective action method 

A second large- expansion is provided by a 2PI effective 
action approach that involves a diagrammatic expansion 
of the double Legendre transform r[?/', C] of the generat- 
ing functional of connected correlation functions (Valageas 
2007a). This is also the generating functional of two-particle 
irreducible (2PI) diagrams, that is those that cannot be 
disconnected by cutting two lines. This approach yields the 
same Schwinger-Dyson equations (|27|) - (p8|) . but the self- 
energy is now given in terms of the nonlinear two-point 
functions C and R. Therefore, this method does not in- 
troduce the auxiliary functions Co and i?o, but the sys- 
tem (|27I) - (P5)) is now a pair of nonlinear coupled equations 
(S oc RC and 11 cx GC at one-loop order). As shown in 
Valageas (2007a), this nonlinearity is the key to the small- 
scale damping of the nonlinear response R (at one-loop or- 
der). At one- loop order, the self-energy has the same struc- 
ture as for the direct steepest-descent method and it reads 
as 

E(a;, y)^AKs(x\ xi,X2)Ks{z; y, Z2)R{xi, z)C{x2,Z2) (31) 
^{x, y)^2Ks{x; xi,X2)Ks{y; yi,y2)C{xi,yi)C{x2,y2) (32) 

Thus, Eqs.dSIl)-® are equal to Eqs.dlHD-dSni) where the 
linear two-point functions are replaced by the nonlinear 
ones. This simple correspondence would not hold at higher 
orders as the self-energies associated with the steepest- 
descent approach would involve additional terms as com- 
pared with the 2PI effective action approach. 

Approximation schemes based on equations of the form 
(HZIl-dMl), (ED)-®, have been apphed to many fields of 
theoretical physics, such as quantum field theory (e.g. Zinn- 
Justin 1989; Berges 2002). They can be derived in differ- 
ent manners and are also called "mode-coupling approxi- 
mations" , in condensed matter and studies of glassy sys- 
tems (Bouchaud et al. 1996), or "direct interaction approx- 
imation" in studies of turbulent flows (Kraichnan 1961). A 
derivation inspired from works on turbulence is presented 
in Taruya & Hiramatsu (2007). An advantage of the deriva- 
tion from the path integral (|24p is its relative simplicity and 
the well-known diagrammatics associated with the expan- 
sion up to high orders. Moreover, it satisfies an exact vari- 
ational principle (the equations above are associated with 
the saddle-point of the 2PI effective action F and it is this 
action which is approximated by truncating an expansion 
over 1/N, see Valageas 2004). Although we do not use this 
property here, it might prove useful at some point. On the 
other hand, the path integral may serve as a basis for 
other approximation schemes. 

2.6. Comparison between various approaches 

The standard Schwinger-Dyson equations l|27|) - ([28|) can also 
be obtained in a diagrammatic fashion from the IPI ef- 
fective action r['0], also called the generating functional of 
proper vertices (see any textbook on quantum field theory) . 
It is the generating functional of one-particle irreducible 
(IPI) diagrams, that is those that cannot be disconnected 
by cutting one fine, and it is also the Legendre transform of 
W = lnZ, where Z is given by Eq. ([M)) . Thus, Scoccimarro 
& Crocce (2006a) describe a diagrammatic derivation of 
Eqs. ((27|) - ((28|) . starting from the equation of motion ([8]) 
written in integral form (as in Ea. (|79p below). 



As noticed above, Eas. p7)l -(|28 p alone are mostly a 
rewriting of the problem, and one still needs to specify a 
method to compute the self-energies S and 11. We have de- 
scribed above two such methods, the direct steepest-descent 
approach, where the self-energy is written in terms of the 
linear two-point functions Cl and Rl, and the 2PI effec- 
tive action approach, where the self-energy is written in 
terms of the nonlinear two-point functions C and R. Of 
course, both methods can also be obtained by diagram- 
matic means. For instance, Scoccimarro & Crocce (2006a) 
present a loopwise expansion in terms of the nonlinear two- 
point functions that coincides with the 2PI effective action 
approach (see also L'vov & Procaccia 1995 for a detailed 
description of such a diagrammatic procedure). 

Next, it is possible to reexpress the self-energy in a more 
compact fashion by introducing "dressed vertices" that go 
beyond the tree-order vertex Kg. That is, one can recognize 
three-point diagrams such as those of Fig. [TO] (but writ- 
ten in terms of nonlinear two-point functions) in parts of 
larger diagrams, and identify them with the expansion of 
properly defined dressed vertices (L'vov & Procaccia 1995). 
However, this leads to another expansion scheme than the 
strict 2PI scheme, that would be more closely related to the 
IPI effective action approach. 

Indeed, looking for an expansion of the self-energies S 
and n in terms of the nonlinear R and C does not uniquely 
specify the practical procedure used to compute correlation 
functions. One can always reorganize the diagrammatic se- 
ries in different fashions that lead to different truncations. 
For instance, as discussed in Sect. 5 of Valageas (2004), 
and recalled above, it is possible to derive the Schwinger- 
Dyson equations ((27)) - (|28p from the IPI effective action 
F [?/>]. Then, the self-energy is given by an exact expres- 
sion that involves the proper vertex F3 which "renormal- 
izes" the bare vertex Ks (its first two terms are given in 
Eqs. ((7D)) - ((75)) below). Next, the proper vertices F, obey a 
hierarchy of equations, that must be truncated at some or- 
der q (just like the familiar BBGKY hierarchy obeyed by 
the correlations Cq can be truncated at some order to build 
an approximation scheme in the weakly nonlinear regime). 
This approach provides another expansion scheme, written 
in terms of nonlinear two-point functions, that also con- 
tains infinite partial resummations as compared with the 
standard perturbative expansion recalled in Sect. [3] below. 

For completeness, let us point out that we always con- 
sider the limit 77/ — > —00 (i.e. the initial conditions are 
set up at a redshift zj —>■ 00) and the response function 
R{xi,X2) is defined as the deviation from the mean of the 
fields at time rji with respect to a perturbation at an arbi- 
trary earlier time r]2 - By contrast, the calculations presented 
in Scoccimarro & Crocce (2006a, b) consider the system de- 
fined by initial conditions set at a finite time r/j (which 
defines the origin of time as rji = 0) and focus on the 
propagator G defined from the restriction of i?(a;i,a;2) to 
r]2 = r]i. However, the formalism of Scoccimarro & Crocce 
(2006a, b) can also be applied to the case 77/^—00 (e.g. 
Sect. [6] below). 

2.7. Zeldovich dynamics 

The gravitational dynamics described in the previous sec- 
tions is rather complex and analytical computations of- 
ten lead to complicated expressions. Therefore, it can be 
interesting to investigate a closely related dynamics that 
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obeys a similar equation of motion but has a simpler struc- 
ture. The Zeldovich dynamics offers such an opportunity 
(Zeldovich 1970; Gurbatov et al. 1989). As described in 
Valageas (2007b), this dynamics can also be defined from 
the hydrodynamical equations (HJ-lHl), where the velocity 
potential is used in place of the gravitational potential. 
Then, the equation of motion again has the form ([8]). The 
quadratic vertex Kg is unchanged whereas the linear oper- 
ator O is changed to 

0{x,x')^(^^ g^\^5D{^-\^')5D{v-v')- (33) 

This yields the same linear growing mode (|16p hence the 
same linear two-point correlation (jlSp . However, the linear 
decaying mode is modified by the change of O; this also 
modifies the linear response as 



Rl{xi,X2) = 5z)(ki - k2) 9{rii - 772) 
lA , / 1 -1 



1 







(34) 



Then, the path-integral formalism of Sect. 12.31 can be ap- 
plied in exactly the same manner and we obtain the same 
action (|25p . where we must only substitute the new linear 
operator O (Valageas 2007b). Therefore, large- iV expan- 
sions again lead to Eqs. ((7f|) - ((5^ . Thus, the Eulerian equa- 
tions of motion have the same structure and any expansion 
scheme can be applied in identical manner to both dynam- 
ics. However, in Lagrangian space the Zeldovich dynamics 
is much simpler as its solution is given by 



X = q-l- i:)+(77)sLo(q) 



(35) 



where q is the Lagrangian coordinate of the particle. That 
is, the trajectories of particles are given by the linear dis- 
placement field Si = _D+Sio- Going back to Eulerian space, 
one can derive from the solution psp an explicit expression 
for the matter density contrast. Starting from the uniform 
matter density p at i ^ 0, the conservation of matter gives, 
before orbit-crossing. 



p(x)dx = pdq whence 1 -I- 5(x) = 



det 



9x 
9q 



(36) 



This also reads from Eg. ([55]) (Schneider & Bartelmann 
1995; Taylor & Hamihon 1996) as 



Si^,v)^ J dq (5i3[x-q-si,(q,?;)] -1. 
In Fourier space we obtain 



dq 

(2^ 



-ik.q 



-ik 



si,(q,i7) _ I 



(37) 



(38) 



whereas the linear displacement field is related to the linear 
density contrast by 



(39) 



From Eg. ([38)1 one can also derive the exact expressions 
of the nonlinear correlation C and response R (Schneider 
& Bartelmann 1995; Taylor & Hamilton 1996; Valageas 
2007b). 



3. Standard perturbation theory 

In this section, we recall the standard perturbation the- 
ory applied to the cosmological gravitational dynamics 
(Fry 1984; Goroff et al. 1986; Bernardeau et al. 2002; 
Scoccimarro 1997). As explained in Sect. 12.71 note that all 
expressions derived below apply as well to the Zeldovich 
dynamics. In this approach, one expands the density and 
velocity fields over powers of the linear growing mode (|16p 
as {tpi = S) 

00 „ 

<5(k, D) ^Y.^" dqi..q„5i3(qi + .. + q„ - k) 

xF„(qi, ..,q„)(5Lo(qi)--'5Lo(q«), (40) 
and = -V.w/Uf) 

M^D) = dqi..q„5D(qi + .. + qn-k) 



x£'„(qi, ..,q„)(5Lo(qi)--<^Lo(qn), 



(41) 



with Fi = El = 1 (in Eqs. (gHll-lIlIl) and in the following q^ 
is not a Lagrangian coordinate but a wavenumber associ- 
ated with a linear growing mode). Here (and in most cases 
below) we use the linear growth factor D = e'' of Eq. ^ as 
the time coordinate. The factorization of the time depen- 
dence as I?" is due to the approximation Om//^ — 1 (and it 
is exact for the Einstein-de Sitter cosmology). Then, by sub- 
stituting the expansions ([40 |) - (|4T|) into the equation of mo- 
tion ([8]) one obtains a recursion relation between (i^„,i?„) 
and (Fn+i, En+i). This enables one to compute the kernels 
Fn and the density and velocity fields up to the required or- 
der. Moreover, one can easily check that the kernels F„ and 
En obey the symmetry F„(-qi, .., -q„) = K(qi,..,q„) 
since (5(x) and 0(x) are real (whence (5(— k) = (5(k)*). They 
also satisfy i^„(qi, .., q„) — for qi + .. + q„ = as 
(5(k = 0) = at all times because of the conservation of 
matter. It is convenient to introduce the symmetric ver- 
tices F^ obtained by summing over all permutations of 
(qi, ..,q„) as 



F,^(qi,..,q„) 



1 



perm. 



^^„(qi,. 



(42) 



5 = 



-0 + 





+ ... 



Fig. 1. The diagrams associated with the standard pertur- 
bative expansion (PD|) of the density field (5(k, rj) over powers 
of the linear growing mode 6l- The filled circles are the ver- 
tices Fn attached to n linear modes 6l shown by the open 
circles. The left dotted legs are attached to the external 
wavenumber k of the nonlinear density (5(k, 77). 

The diagrammatic expansion associated with (|40p is 
shown in Fig. [TJ The diagrams are very simple to draw 
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but the complexity is hidden in the vertices that are 
obtained by recursion from lower-order vertices. Note that 
one needs to compute a new vertex at each order of the 
expansion. 

Then, from the expansion (j40p . one obtains for any cor- 
relation function a perturbative expansion over powers of 
the linear power spectrum Pl by using Wick's theorem. In 
particular, the two-point density correlation Cs2 reads (up 
to order whence up to order P|) as 



C, 



52 



(66), = -I- 

+(<5(2)<5(2^) + ..., 



(^(1)^(3)) 



(43) 



where J*^"^ oc F^S^ is the term of order n in the expansion 
(|40|) . The diagrammatic expansion associated with Eq. ((43l) 
is displayed in Fig. [21 The order of the diagrams over 
scales with the number of loops £ as P^'^^. Up to order P|, 
we can write the equal-time nonlinear power spectrum P{k) 
as 



P{k) = P^"'%k)+P^^°°P{k), 
with 

ptioc _ p(a) _ P]^{k) pll°°P = pO>) 



(44) 



(45) 



Thus, the tree diagram (a) of Fig. [2] is merely equal to the 
linear power spectrum P^ whereas the one-loop contribu- 
tions pf'') and of diagrams (b) and (c) are given by 



pW(A) = 6Pi(fc) / dqFi(g)P|(q,-q,k), 



F(=)(fc) 



2 / dqPi(g)Pi(|k-q|)P|(q,k-q)2 



(46) 
(47) 



We give in Fig. [2] the overall multiplicity factor of each di- 
agram. Thus, the factor 6 for diagram (b) in Fig. [2] reads 
as 6 = 2 X 3, where the factor 2 comes from Ea. (|15)l (as- 
sociated with the two contributions (J^^^J^^^) + ((J^^^J^^^)) 
whereas the factor 3 comes from the three possible Gaussian 
pairings which can be built for each term (J^'^^tJ^^^) (Wick's 
theorem). In this article, we focus on equal-time statistics, 
such as the usual matter power spectrum P(fc), and each 
factor of Pl gives rise to a time-dependent factor D^, where 
D is the linear growth factor at the time of interest, but the 
analysis applies as well to different-time correlations. 



(a) 



+ 6 



(b) 



+ 2 



(c) 

o - 



Fig. 2. The diagrams associated with the standard pertur- 
bative expansion (j43p of the two-point density correlation 
Cs2 over powers of the linear power spectrum P^. The big 
dots are the vertices F^ and the solid lines are the linear 
power spectrum Pl- We also give the overall multiplicity 
factor of each diagram, up to one-loop order (i.e. P| for 

CS2)- 

On the other hand, the three-point density correlation 
Cs3 reads (up to order P|) as 

C53 = ((5(2)<5(i)5(i))-f 2perm. + (^(^)<5(i)(5(i)) + 2perm. 

+ (<5(3)j(2)^(i)) + 5pcrm. + (<5(2)5(2)^(2)) + ^ (43) 



(a) , ■ (b) • (c) 
C.3= 6 <]+36 0+36 

(d) K (e) 

+ 36 +8 





+. 



Fig. 3. The diagrams associated with the standard pertur- 
bative expansion (pS)) of the three-point density correlation 
Cs3 over powers of the linear power spectrum Pl (shown 
up to one-loop order, i.e. P| for C^s). The symbols are as 
in Fig. H 



where each term "n perm." stands for n terms obtained by 
permutations over the previous term, such as 

(^(2)^(1)^(1))+ 2 perm. = (5(2)^(1)^(1)^^^^(1)^(2)^(1)^ 

+ (5(1)^(1)^(2)). (49) 

The associated diagrams are shown in Fig. [31 For the equal- 
time bispectrum B defined by 

C'53(ki,k2,k3;?7) = fo(ki + ka + k2)P(fci, ^2, /cs; ??), (50) 
we obtain at tree-order the well-known result 
P"'=°(fci,fc2,fc3)- 2Pi(fc2)-PL(fc3)-F|(k2,k3) + 2perm. (51) 
and at one-loop order: 

P"°°P(fci, A:2, fc3) = S^*") + P(") + P^'') + B^'\ (52) 
with (Scoccimarro 1997): 

P(^) = l2PL{k2)PL{k3) y"dqPL(g)P4^(-q,q,k2,k3) 

-|-2perm., (53) 

P(^) = 6PL(fc2)-PL(fc3)-F2(k2,k3) J d(lPL{q) 

xP3(q, -q,k2) + 5perm., (54) 

P^"^) = 6PL(fci) ydqidq2 ^^(qi + q2 - k2)PL(<Zi) 

xPL(g2)i^|(qi>q2)-F'3(qi>q2,ki) + Sperm., (55) 

P^"' = ^ J dqidq2dq3 ^^(qi -|- q2 - ki) 

xSoicii + qs + k2)PL{qi)PL{q2)PL{q3) 

X P| (qi , q2 )P| (qi , qs)^^! (-q2, qs) ■ (56) 

The derivation of higher-order correlations proceeds in 
exactly the same manner, by writing Csp = {Si..Sp), ex- 
panding each factor Sj through Ea. (|40p . and performing 
the Gaussian average with Wick's theorem. 
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4. Expansion over cubic interaction term 



R, (x,,x,) 



In this section, we recall how to recover the results of the 
standard perturbation theory from the path integral (j24p . 
As described in Valageas (2004, 2007a), the usual pertur- 
bative expansion presented in Sect. [3] over powers of the 
linear growing mode Sl can also be seen as an expansion 
over powers of the interaction vertex Ks of the equation of 
motion ([8]). (Formally one may add a coupling constant g as 
Kg —> gKg and expand over powers of g.) Indeed, by sub- 
stituting the expansions (pn)l - (liT|) into the equation of mo- 
tion ([8]) to derive the recursion relation between (F„,i?„) 
and (F„+i, i?„+i), we can see that Fn oc JsT""^ whence 
J*^"-* oc K^~^ (by which we actually mean (5^"-' oc g""^). 
Therefore, the expansion over powers of 6l, or over powers 
of Pl for averaged quantities such as correlation functions, 
is identical to the expansion over powers of Kg. Then, this 
latter expansion can be directly obtained from the path in- 
tegral by expanding over the cubic part X.Ks-tpip of 
the action S[ip,X] of Eq. ljSS]) . This equivalence applies to 
all higher-order correlation functions. 

Thus, we obtain for the two-point function 



'[d^][dA]^(a;i)V'(x2) [l + XKs^^ 
= {^P{xi)^P{x2) (l + XK,^j^ + 



2 , (Ai^.^^)2 



2 

+ -))o 



-So 



(57) 



where 5*0 — X.O.iIj — ^A.A/.A is the quadratic part of the 
action S'['0, A] and (..)o is the average with respect to this 
Gaussian action 5o. Then, one can use Wick's theorem to 
compute the Gaussian path integral defined by 5*0. Since 
the action is equal to the action 5* with = 0, it 
actually corresponds to the linear dynamics. Thus we have 
(as can also be explicitly checked): 

(V')o = (A)o = 0, (^(a:i)^(x2))o = Cl{x^,X2). (58) 
(^(a:i)A(a;2))o =i?L(.Ti,X2), (A(.ti)A(x2))o = 0. (59) 

In terms of diagrams, we note these linear propagators as 
in Fig. m We add an arrow to the lines associated with 
the propagator Kh — {'<jj{xi)X{x2))o to mark the arrow of 
time due to causality {Di > D2), which is enforced by the 
Heaviside factor of Eg. p^ . Thus, the response Rl runs 
from the response-field A towards the physical field ip. We 
could have added an out-going arrow at both ends of Cl 
(both ends are associated with fields tp) but for simplicity 
we put no arrow at all. Since each vertex Kg is associated 
with one response-field A and two physical fields '0, through 
XKsipip as in the action (PS)) . each vertex Kg shown by a 
big dot in Fig. [S] and hereafter must be connected to one 
outgoing line (hence a response Rl) and to two incoming 
lines (which can be Rl or Cl with the understanding that 
Cl can be read with two outgoing lines, that is, one at each 
end). 

Then, from Eq. (|57p the two-point correlation reads up 
to one-loop order as 



^trcc _|_ ^lloop 



with Cr'' = Cl, 



(60) 



and 



Fig. 4. The diagrammatic symbols used for the two-point 
functions Cl = (V''0)o and Rl = {'4'X)o- The third possi- 
ble two-point function vanishes: (AA)o = 0. The propaga- 
tor Rl = {ip {xi)X{x2))o is distinguished from Cl by the 
addition of an arrow which marks the direction of propa- 
gation, from the response-field A(x2) towards the physical 
field ip{xi), that is, from time D2 to time Di > D2 as en- 
forced by the Heaviside factor of Eq. (fT9|) . (This is due to 
causality.) 



8RLiKgRLCLKg)CL ^ 

+mLCLKg{RLKgCL) 



2RL{KgCLCLKg)RL 



(61) 



The diagrams associated with Eqs. ([50)) - ([CT|) are shown in 
Fig. [5] In Eq. ([6T|l and in Fig. [5] we have not shown five 
additional terms and diagrams that vanish because they 
contain a two-point function (AA)o or a closed loop over re- 
sponse functions Rl- Indeed, the former vanishes because of 
Eq. ([55)) (and more generally Eq. (pS|) ) whereas closed loops 
over Rl vanish because of the Heaviside factors 6* (771 — 772) 
associated with Rl as in Eq. p9)) . Moreover, it happens 
that diagram (d) of Fig. [5] also vanishes because the right 
vertex Kg attached to the closed loop over Cl is of the 
form Kg{Q; q, — q) (where we only write the dependence on 
wavenumbers) , which is equal to zero from the explicit ex- 
pressions ((12)- (HH) of the vertices 7". Then, one can check 
that the diagram (b) gives back diagram (b) of Fig. [2] 
and Eq.(f46|) whereas diagram (c) gives back diagram (c) 
of Fig. [5] and Eq. ([T7)) . Thus, we indeed recover the stan- 
dard perturbative results. Nevertheless, these two equiv- 
alent expansions have different diagrammatic structures. 
Thus, the expansion of Fig. [51 derived from the path inte- 
gral ((24)) . only involves the cubic vertex Kg but two prop- 
agators Cl and Rl (with four indices) because of the two 
two-component fields ip and A. By contrast, the standard 
expansion of Fig. [2]involves new vertices of increasing or- 
der as we include higher-order terms but only one two-point 
function: the linear density power spectrum Pl- However, 
all kernels Fn can be written in terms of the building blocks 
Rl and Kg, as in Fig. [13] and Eqs.(l79l)-(l80l) below, which 
are also the basic blocks of the expansion in Fig. [5] 



(a) 



+ 2 




+ 8 



(c) 




(b) 



(d) 



C 



Hoop 



-(V(xi)7/.(x2)(Aif,^V)2)o 



Fig. 5. The diagrams associated with the expansion ([57| 
over powers of Kg for the two-point correlation C2. The 
big dots are the three- leg vertex Kg . The solid lines are the 
linear propagators Cl and Rl as in Fig. [H Note that the 
symbols are different from those used in Figs. [TJH for the 
standard perturbative expansions. Again we only display 
the diagrams obtained up to one-loop order (i.e. P| for 
C2). 
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Higher-order correlation functions can be obtained as 
in Ea. (|57p by expanding over Ks- This gives for the three- 
point correlation: 

{^{xi)i;{x2mx2) (\K,^? + + . Vq. (62) 



6 

Thus, we obtain up to 1-loop order 
CT'= = i?Li^sC^LCL-t-5perm., 
and 

Cl^°°^ ^ (24 -f 48 -f 48){K^RIcI)KsRlCl 

3r<3 



(63) 



ct) 

A8)RlCI{kIrICl). (64) 



+2A(KIrICl)KsCI + 8i?i(X 

+AmlCL{KlRLCl) + 

We show the diagrams associated with Eas. ([63| - (|64p in 
Fig. [HI We can note that the first four one-loop diagrams, 

(b) , (c), (d), and (e), are 1-particle reducible: they can be 
disconnected by cutting the two-point function which con- 
nects the "bubble" to the rightmost vertex Kg- The com- 
parison with Fig. [5] shows that they correspond to the per- 
turbative corrections of the two-point functions. Thus, the 
first three one-loop diagrams can be obtained from the tree 
diagram (a) by inserting the one-loop corrections (b) and 

(c) of Fig.[5]into the two-point correlations Cl whereas the 
fourth one can be obtained by inserting the one-loop correc- 
tion of the response R. The last four one- loop diagrams (f), 
(g), (h), and (i), of Fig. [6l correspond to perturbative cor- 
rections to the vertex Kg itself. Note that whereas the last 
two diagrams respect the structure of the bare vertex Kg 
(attached to one response field A and two fields tp) the other 
two diagrams bring up vertices attached to three fields A 
or to two fields A and one field ip- Thus all possible con- 
figurations are generated beyond tree-level for the "renor- 
malized" vertex. As for the two-point function diagrams 
shown in Fig. [5l we have only displayed in Fig. [5] the non- 
vanishing diagrams. For instance, the renormalized vertex 
which would arise from a closed triangular loop over three 
responses Rl vanishes because of the Heaviside factors as- 
sociated with the response Rl (that express causality). 

Again, one can check that Eas. ([55|) - ([M|) give back 
the results obtained from standard perturbation theory in 
Eq.lgll) and Fig. H Thus, from the expression (HH) of the 
bare vertex Kg, the tree-order diagram reads as 

rvi 

Cf'='=(a:i,a;2,a;3) = (5i3(ki +k2 +k3)/ d-q'iRLhi'^[ki;T]i,ii^) 

J —OO 

X CL^2^'^ (fc2 ; »72 , v'l)Cu3t'^ (h^VS, v'lllt'^ {^2 , k3) 

-I- 5 perm. (65) 

Substituting the expressions of the linear two-point func- 
tions in the right hand side of Eq. (|65p and performing the 
integral over the time rj'^ we recover Eg. ljSip for ii =12 — 

»3 - 1. 

Higher-order correlations Cp are computed in the same 
fashion, by inserting the factor {'ipi..^p) in front of the ex- 
ponential in Ea. ([M)) and expanding over the cubic part of 
the action. The same technique also provides higher-order 
response functions, such as 



{■)p{xi)->p{x2)) 



6C{xi,X2) 

sax3) 




(b) 



+ 48 ^ 



+ 24 ^ 



+ 48 



+ 8 



+ 48 




+ 24 ^ 




(66) 



+ 48 



Fig. 6. The diagrams associated with the expansion ([62|) 
over powers of Kg for the three-point correlation C3. The 
symbols are as in Fig.O Again we only display the diagrams 
obtained up to one-loop order (i.e. P| for C3). 

Thus, i?0^A measures the response of the two-point corre- 
lation C to an infinitesimal external perturbation ( added 
to the right hand side of the equation of motion ([H]) . It is 
also related to the cross-correlation of the nonlinear field tp 
with the initial condition ijjj at a finite time 77/ as 

{ipixi)ip{x2)'ipiix3)) ^ R^^i,x{xi,X2;x,rii) X Gi{x,X3), (67) 

where the product does not involve any integration over 
time and we note x = (k, i), see Eq.(49) of Valageas (2007b) 
for details in the case of the two-point response function. 

Note that throughout this article we take the limit 
rji —00 (i.e. Dj — > 0). This applies both to the standard 
perturbation theory of Sect. [3] (and the expansion (|79p be- 
low), which allows us to write the solution of the equations 
of motion as an expansion over the linear growing mode, 
and to the methods derived from the path integral (|24p . It is 
possible to apply the path-integral formalism while keeping 
■qi finite, see Valageas (2007b), but computations simplify 
in the limit 77/ — > — cx). In a similar fashion, keeping 77/ finite 
would add decaying terms to the standard expansion (|40p . 



5. Large- methods 

We now describe in this section how to obtain the three- 
point density and velocity correlations from the large-A^ 
methods introduced in Sect. 12.31 From the path- integral 
formalism (|24p - ([^5|) the three-point correlation C3 can be 
written as (Zinn- Justin 1989; Valageas 2007a) 

C3(a;i,a:2,a;3) = -C{xi, x[)C{x2, X2)C{x3, x'3)r3{x[, x'2, x'^) 

(68) 
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where T3 is the symmetric dressed three-point vertex (the 
so-called three-point proper vertex of the IPl effective ac- 
tion r['0]). Here we put on the same footing the fields "0 and 
A so that two-point functions C correspond to all three pos- 
sibilities (tpip) = C, {tpX) = R, and (AA) = 0, and similarly 
for the three-point function C3. Then, the proper vertex Fg 
can be obtained from a series of three- leg diagrams, up to 
the required order over 1/N. 

In this article we restrict ourselves to the three-point 
correlation C3{xi, X2, x^,) defined as 

C3(a:i,a:2,a;3) = {il:{xi)ip{x2)il^{x3,)) . (69) 

It describes both density and velocity correlations, as well 
as their cross-correlations. We do not study in this paper 
three-point functions such as {ipipX) of Ea. ((66)) . associated 
with the functional derivative of the two-point correlation 
with respect to an external noise (and in particular with 
respect to the initial conditions). Of course, higher-order 
correlation and response functions can be obtained in a sim- 
ilar manner, from the two-point function C and R and the 
higher-order proper vertices Tp. However, since for practi- 
cal purposes higher-order statistics are increasingly noisy 
we do not go beyond three-point functions in this paper. 

Note that the structure of the large-A^ expansions 
is somewhat different from the expansions described in 
Sects. [3]|31 Indeed, contrary to the latter cases, the contri- 
bution associated with a given diagram is not fixed but de- 
pends on the order up to which the expansion is performed. 
This arises from the fact that each diagram depends on the 
nonlinear two-point functions R and C (more precisely the 
approximated nonlinear two-point functions computed at 
this order). Thus, the diagrammatic series actually gives 
an implicit equation for R and C, as in Eas. ([?71) - (P5)) where 
the right hand side depends on R and C. 

Then, the order of these expansions is set by the or- 
der of the diagrams kept in the self-energies E and H. At 
tree-order, we merely have E = and H = 0, at one-loop 
order we have either Eqs.(l29l)-(l30l) or Eqs.(IST|)-(l32l), and 
at higher orders we would need to include higher-order dia- 
grams (that would no longer look identical for the steepest- 
descent and 2PI methods). 

5.1. Direct steepest-descent method 

We first consider in this section the direct steepest-descent 
method of Sect. 12.41 At tree-order, the two-point functions 
are equal to the linear two-point functions (Valageas 2004, 
2007a) whereas the proper vertex Fa is equal to the bare 
vertex of the action S[tp,X] which yields from Eq. ipS)) 

F^''°'=(a::i,a;2,a;3) = -Ks{xi;x2,X3) + 5 perm., (70) 

with a multiplicity factor 6 = 3! which corresponds to all 
permutations of the triplet {xi,X2,X3). Therefore, at tree- 
order we recover exactly the tree-diagrams (a) of Figs. [5]- 
[SI obtained by expanding over the nonlinear coupling K^, 
hence the lowest-order results of standard perturbation the- 
ory 

At one-loop order, the two-point functions are obtained 
from the system (HTJ-dlHl) together with Eqs.(I211)-(ISni) for 
the self-energy. The diagrammatic expression of the self- 
energy is shown at one-loop order in Fig. [71 Then, the so- 
lution of Eq. (1^51) can be written as 




Fig. 7. The one-loop diagrams for the self-energies E and 
n obtained from the direct steepest-descent expansion. 




Fig. 8. The diagrams obtained at one-loop order for the re- 
sponse R from the direct steepest-descent expansion. The 
first equality corresponds to the implicit equation (|7ip . 
where we note by double lines the nonlinear response R. 
The second equality shows the infinite series of bubble di- 
agrams resummed by this large- method, as in Eg. ([7^. 

= Rl + Rl-^-Rl + Rl-^-Rl-^-Rl + - (72) 

In the first line we used the second Eq. p5|) to obtain the 
integral form of Eq. ([28|) and in the second line we wrote 
the solution as a series over powers of E. Equation (|72|) ex- 
plicitly shows that at "one-loop order" the steepest-descent 
method has performed the resummation of the infinite se- 
ries of bubble diagrams displayed in Fig. [H 

Let us stress here that for the large- A^ methods, the loop 
order of the approximation refers to the truncation order of 
the infinite series of diagrams obtained for the self-energies, 
and not to the truncation order of the correlation and re- 
sponse functions. Indeed, any finite order in terms of the 
self-energy diagrams automatically leads to an infinite par- 
tial resummation over all orders for the diagrams associated 
with the correlation and response functions. However, these 
resummations are only complete up to the same order as for 
the self-energy (we miss some of the higher-order diagrams) . 
Thus, at one-loop order, taking into account the only two 
diagrams of Fig.[7]for the self-energies actually gives rise to 
the infinite series shown in Figs. [MS for the physical two- 
point functions R and C, which include contributions that 
contain an arbitrary number of loops. However, these par- 
tial resummations are complete only up to one-loop, that is 
they miss some two- loop and higher-order diagrams, such 
as the one of Fig. [14] for C . 

Next, the solution of Ea. (|27)) can be written as (Valageas 
2004, 2007a) 

C = i? X Cl(?7/) X i?^ + R.n.R^, (73) 

where the first product does not contain any integration 
over time, 

R X Cl(77/) X R^ ^ (5£i(ki +k2)^i?iii'^(fci;?7i,7y/) 

xCLi[i'^ {ki;vi, Vi)Ri2i'2 iki;V2, Vi), (74) 

and we let the initial time go to the infinite past 77/ — > —00. 
Moreover, we can note that in the linear regime we have 



R — i?L + R^.Y^.R 



(71) 



CL = RLXCLim)xRl, 



(75) 
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Fig. 9. The diagrams obtained at one-loop order for the 
two-point correlation C from the direct steepest-descent ex- 
pansion. The first equahty corresponds to Eq. ()73|) . where 
we note by double lines the nonlinear response R. The sec- 
ond equality shows the two infinite series of bubble dia- 
grams resummed by this large- iV method, in terms of the 
linear functions Cl and Rl- 

as can be checked from the explicit expressions of Cl and 
Rl. From Eqs.(I711)-(17SD and Fig. [5] we see that the infi- 
nite series of diagrams resummed by the steepest-descent 
method at one-loop order are the two series of bubble di- 
agrams shown in Fig. [5] We can check that we recover 
the nonzero diagrams obtained at one-loop order in Fig. 
within the perturbative expansion over the cubic interac- 
tion Kg- In addition, the implicit equations (|77|) - (P5|) have 
allowed us to resum two infinite series of higher-order dia- 
grams. Thus, we recover the fact that the large-iV expan- 
sions and the standard perturbative expansions agree up 
to the truncation order (here one-loop) and only differ by 
higher-order terms (Valageas 2004, 2007a). 

Note that we do not recover the diagram (d) of Fig. O 
that was equal to zero. In fact, the property Ks{0; q, — q) = 
that sets this diagram to zero has already been used twice 
to derive the Schwinger-Dyson equations ((?7|) - ([5(I)) : i) to ob- 
tain the path integral where we used that the Jacobian 
is an irrelevant constant, and ii) to obtain Eqs. (P7|) - (P0)l by 
expanding over 1/N where we used that (ip) = (which 
also comes from Ks{0;q, ~q) = 0). Then, the explicit use 
of these simplifications has already removed from the large- 
N expansion diagrams such as the diagram (d) of Fig. [5] 

For the three-point correlation C3 we need the one-loop 
contribution to the three-point vertex Fa which reads as 
(e.g. Zinn- Justin 1989; Valageas 2004) 

J.1100P ^ ^Q^KsCoKsCoKsCo, (76) 

where the two-point functions Co form a loop which joins 
the three bare vertices Ks- Here Co is the auxiliary two- 
point function introduced within the direct steepest-descent 
scheme as in Sect. 12. 4[ which happens to be equal to the 
linear two-point function Cl- Then, from Eg. ([68)) the "tree- 
diagram" contribution reads as (using Eq. fTDl) ) 

Cl'""" = 6CCC.Ks, (77) 

whereas the "one- loop diagram" reads as (using Eg. ([75])) 

C3"°°P = 216CCC.iKsCLKsCLKsCL). (78) 

We show in Fig. [TO] the diagrams associated with Eqs. ([T7|) - 
([78]) that give the three-point correlation at one-loop order 
within the steepest-descent expansion. In both equations 




Fig. 10. The diagrams obtained at one-loop order for the 
three-point correlation C3 from the direct steepest-descent 
expansion. The double lines are again the nonlinear re- 
sponse R and the nonlinear correlation C. We use the same 
labels as in Fig.[B| The four diagrams (b) to (c) of Fig.[n]do 
not appear in this expansion scheme since their contribu- 
tion is included in diagram (a) through the use of nonlinear 
two-point functions instead of linear ones in the external 
legs. 

0'' 




Fig. 11. An infinite series of diagrams resummed by the 
"tree-diagram" of Fig.[Tn]of the steepest-descent expansion 
at one-loop order. It is obtained from the diagrammatic 
expansions of the nonlinear R and C shown in Figs. [8][9l 
There are two more series (not shown) obtained by taking 
into account the second diagram of Fig. [Hi 

()77p - ()78p we must use the nonlinear two-point functions C 
computed at the same order within the steepest-descent 
method (here one-loop). That is, from the one-loop dia- 
grams of Fig. [7] for the self-energy, which lead to the in- 
finite loop resummations for two-point functions shown in 
Figs.[S][n] (Let us recall that the "loop-order" refers to the 
diagrams kept for the self-energy, but the two-point func- 
tions R and C contain resummations over diagrams that 
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contain an arbitrary number of loops over the linear prop- 
agators.) 

The three nonlinear two-point functions C that appear 
in Eas. (|77|) - ([75)) are shown in Fig. [TU] by the double lines 
associated with the three external legs, as in Eg. ((68)) . Thus, 
the contribution ([77|) . shown by diagram (a) in Fig. [TUl is 
no longer equal to the standard tree-diagram (a) of Fig. [5] 
associated with Eqs. ljGS]) and ([65]) . It includes three infinite 
series of additional diagrams, such as the series shown in 
Fig. 1111 which are obtained by replacing the nonlinear two- 
point functions R and C by their diagrammatic expressions 
shown in Figs. [Sjni In particular, these series contain the 
first four one-loop diagrams (b), (c), (d), and (e), of the 
expansion over Kg shown in Fig. [6l The other diagrams 
of Fig. [TU] also contain such infinite series once written in 
terms of the linear two-point functions Rl and Cl- Thus, 
the steepest-descent method has allowed us to resum in an 
automatic manner some partial infinite series of diagrams 
for all many-body correlation functions. 

5.2. 2 PI effective action method 




Fig. 12. An infinite series of diagrams included in the two- 
point functions given by the 2PI effective action expansion 
at one-loop order. There are additional infinite series of 
diagrams as compared with the steepest-descent method 
thanks to the "renormalization" of the self-energy terms. 

As recalled in Sect. 12.51 the 2PI effective action method 
also gives the Schwinger-Dyson equations (l27 |) -([28 l) . In par- 
ticular, Eqs. (|7T|l - ([73)) still apply but they are now nonlin- 
ear as the self-energy terms now depend on the nonlinear 
two-point functions as in Eqs. ([3T|) -([32 | . At tree-order we 
again recover the results of the standard perturbative ex- 
pansion. At one loop-order the self-energy is given by the 
same diagrams as those of Fig. [7] except that the internal 
two-point functions are the nonlinear ones (i.e. single lines 
must be replaced by double lines). Then, it is clear that 
the two-point functions obtained from Eqs. ([72|l -([73 l) con- 
tain the same infinite series of diagrams as those derived in 
the steepest-descent method, shown in Figs. [5][S1 but also 
additional series due to the "renormalization" of the self- 
energy terms. For illustration we display in Fig.[T2]a typical 
diagram included in the 2PI effective action rcsummation 
for the two-point functions R and C. 

Next, the three-point proper vertex T3 at one loop-order 
can also be written as Eq. (j76p . where the linear two-point 
functions Cq must be replaced by the nonlinear ones C (i.e. 
all single lines including the internal ones must be replaced 
by double lines) . Expanding back over the linear two-point 
functions Rl and Cl we recover again all diagrams ob- 
tained within the steepest-descent method as well as addi- 
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tional series which include terms such as those associated 
with Fig. [HI 

6. Random advection in the large-A; limit 




Fig. 13. The expansion of the nonlinear field tp over the 
linear growing mode tp l from Eq. (j79[) , up to order i/j'^ . The 
filled circles are the vertex Ks of Eq. ([80| , whereas the white 
circles are the linear growing mode ipL- The numbers are 
the multiplicity factor associated with each diagram. 

As recalled in Sect. [31 the standard perturbation theory 
associated with Ea. (|40)) writes the nonlinear density field 
5{\s.) as a series over powers of the linear density contrast 
6l- The kernels Fn of Ea. (j40p are computed by recursion 
from the equation of motion ([8]) so that each order involves 
a new vertex F„. However, it is clear that the latter can also 
be written in terms of the vertex Kg that appears in Eq. ([5]) . 
Thus, the equation of motion ((S]) for the two-component 
field ip can be solved through the expansion over the lin- 
ear growing mode ipL oi Eq. (|16p as (Crocce & Scoccimarro 
2006a,b; Valageas 2007b) 

^ = + kg^l + 2kl^:l + bkli^i + ... , (79) 

where we introduced the integral vertex Kg defined from 
Ks and the linear response i?^ by 

ks = Rl.Ks (80) 

The diagrams associated with the series (|751) are shown 
in Fig. [T21 see Crocce & Scoccimarro (2006a). Of course, 
Eq. ([79)) and Fig.[T3|are equivalent to the standard pertur- 
bation series of Eq. ljlO)) and Fig. [1] Then, with the aim of 
computing the response function i?, Crocce & Scoccimarro 
(2006a) noticed that it is possible to perform a partial re- 
summation of the diagrams shown in Fig. [T3| in a high-A; 
limit. First, one only keeps the diagrams such as the last 
one in Fig. [131 where one can define a "principal path" con- 
necting a linear mode t/^l to the left root so that all other 
linear modes ipL are directly connected to this path. (This 
is not the case for the other diagram of order ?/;|^ in Fig.[T3l) 
Second, in a high-fc limit one may approximate the vertex 
Ks by 

fc>/ci : Ks{x]Xi,X2) ^ (5_D(k2 - k)5£)(?7i - vi)5D{r]2 - v) 

k ki 

'5^,^2 '5*1,2 • (81) 
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This simplifies the series (j79p which can be resummed to 
give the response function (Greece & Scoccimarro 2006a, b) 



i?(a;i, a;2) = RLixi,X2)e' 



(82) 



where is the variance of the one-dimensional linear 
displacement field slo (equal to the variance of the one- 
dimensional linear velocity dispersion up to a normalization 
factor) given by 



An 

T 



dfcPio(fc)- 



(83) 



As seen in Valageas (2007b), keeping only the "principal 
path" diagrams amounts to approximate the equation of 
motion ([8]) by the linearized equation 

0.iP^2K,^L^, or 4> = 4>l + 2K,4>l4'- (84) 

Then, using the approximation (j8ip . the solution of Eq. 
reads as (Valageas 2007b) 



(85) 



The key simplification comes from the approximation 
(5D(ki + k2 - k) ~ (5£,(k2 - k) in Eq.(l8l]). This makes aU 
the intermediate legs in the simplified diagrammatic series 
factorize as powers of J dq ^ (5i(q, 77') which allows the 
resummation ([55]) . Then, it is easy to derive Eg. ([5^ from 
Eq. (l85l) . as well as 



(86) 



As discussed in Valageas (2007b), this procedure has ac- 
tually replaced the gravitational dynamics by the simple 
effective dynamics 



^(X,Z?) =5L(x-SL(xLag =0,D),D). 



(87) 



That is, the linear density field is uniformly advected by the 
linear displacement SL(xLag — 0,D) associated with the 
Lagrangian position XLag = 0. Then, the average over this 
random Gaussian displacement leads to the apparent loss of 
memory seen in Eas. ([5^ . (|86p . through the Gaussian decay 
^ As discussed in Valageas (2007b), it happens 
that for the Zeldovich dynamics Eg. ([5^ actually gives the 
exact nonlinear response whereas Eg.(|86p gives the correct 
Gaussian decay for different times Di ^ D2 but does not 
capture the equal-time behavior of the nonlinear correlation 
C. 



+ .... 



Fig. 14. A diagram contained in the high-fc resummation 
associated with Eq. ([M)) for the two-point correlation C. 




As compared with the large- resummation of Sect. 15. ll 
which gives the bubble diagrams of Fig. [9] for the two-point 
correlation C, the resummation associated with Eq. (|M|) 
(i.e. keeping the principal path diagrams for ij:) contains ad- 
ditional diagrams such as the one shown in Fig.[T4l Indeed, 



when we join the relevant diagrams of Fig. [T3] to compute 
C(xi,X2) = {4'{xi)ip{x2)), the Gaussian average draws all 
possible pairs from all secondary legs ipL{(lj) that connect 
to the principal path (including those obtained by match- 
ing fields '0L(qi) that are associated with both nonlinear 
fields 'ip{xi) and ip{x2))- By contrast, the bubble diagrams 
of Fig.[5]only involve the Gaussian pairings performed in se- 
quential order as we move along the principal path (but this 
method does not require the high-fc approximation (IST]) '). 

Nevertheless, at one-loop order we can substitute the 
two-point functions ([5^ . ([55]) . obtained in this fashion into 
the diagrams of Fig. [TO] to derive an approximation for the 
three-point function that is correct up to this order and 
contains infinite partial resummations. We can obtain two 
approximations for C3 in this manner. We can either use 
Egs. (|82| . (|86|) . into the external lines only of the one-loop 
diagrams of Fig. [TOl in the spirit of the steepest-descent 
scheme of Sect. 15. 1[ or into both the external and inter- 
nal lines in the spirit of the 2PI effective action method 
of Sect. 15.21 As for previous schemes, these two procedures 
may be applied to all higher-order correlation and response 
fimctions. 



7. Zeldovich dynamics 

We investigate in this section the predictions of the various 
expansion schemes described in the previous sections for the 
simpler case of the Zeldovich dynamics (Zeldovich 1970), 
where the trajectories of particles follow the linear displace- 
ment field as in Eq. ([55]) . We focus on the three-point func- 
tion C3, and more specifically on the matter density bis- 
pectrum defined in Eq. (|50p . since the two-point functions 
C and R have already been studied in Valageas (2007b). 

7.1. Standard perturbation theory 

As is well-known, from the exact solution ([55]) of the 
Zeldovich dynamics the kernels Fn introduced in the stan- 
dard perturbative expansion (|40p can be directly derived at 
any order (Grinstein & Wise 1987). For instance, expanding 
the exponential in Eq. (|55|) gives 



^(k)-|:^/(0i^-'^'^Hk.s.)". 



Using Eq. ((39)) we recover Eg. (|40)l with the symmetric ker- 
nels 



■F;'(qi, --^qn) 



1 k.qi k.q„ 



9? 



In 



(89) 



Then, using Eqs. (|15)l - ([17)) the power spectrum reads up to 
one-loop order as 



P(fc) = PL(fc) - PL(fc) J dq^^^Pi(q) + y"dqidq2 

(k.qi)2(k.q2)2 



^Soiqi + q2 - k)PLiqi)PL{q2 



Hqt 



(90) 



Egs. (|5T|) - (|56)) read for the bispectrum up to one-loop order 



P"-°°(fcl,fc2,fc3)=PL(fc2)PL(fc3) 



;ki.k2)(ki.k3 

"-2 "-a 



-2perm.(91) 
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and, 



^K2K^ J 



(ki.q) = 



+2 perm. 



(92) 



-PL{k2)PL{k3) ^7:212 /dqPL(g) 4 



2 A^2 ^2 



-5 perm., 



(93) 



B^"^ = -PL(ki)^^ I dqidq2 Soiqi + q2 - ka) 
(k2.qi)(k2.q2)(k3.qi)(k3.q2) 



xPl(<Zi)Pl('?2)- 
-1-5 perm.. 



4 4 



(94) 



B^'^) = -ydqidq2dq3 ^^(qi + q2 - ki) 

X Soim + qs + k2)PL(gi)i'L('?2)^'L(g3) 

(ki.qi)(ki.q2)(k2.qi)(k2.q3)(k3.q2)(k3.q3) 



444 

mm 



(95) 



Moreover, we can gather diagrams (b) and (c), which have 
the same form, as 

5W+5(^) = _p,(fc^)p^(fc3)0ii:^K^ /dqP,(,) 
^(k,.q)^ + (k2.q)- + (k3.q)-^^^^^^^^ (96) 



We can check that for a CDM-like linear power spectrum, 
with Phik) oc fc for — > and Phik) oc for k ^ 00, all 
contributions to the power spectrum and to the bispectrum 
converge. More precisely, if we define the local spectral in- 
dex n as 



7.2. Expansion over cubic interaction 

We now turn to the method presented in Sect.[4l where we 
expand the path integral over the interaction vertex Kg. 
As explained in Sect. IH this method exactly recovers the 
results of the standard perturbation theory of Sect. 17.11 

First, as discussed in Sect. [H for the nonlinear power 
spectrum the diagrams (b) and (c) of Fig. [5] give back the 
contributions of diagram (b) and (c) of Fig. [2] associated 
with the standard perturbative expansion. Note that one 
can see at once that there is no mixing of diagrams as 
one goes from one expansion to the other one since it is 
clear that both diagrams (b) contain a fixed factor Phik) 
whereas both diagrams (c) contain a double integral over 
PL{qi)PL{q2) (hence there is a one-to-one correspondence). 
In any case, a simple explicit computation of the diagrams 
of Fig. [5] gives indeed Eg. ([50]) for the Zeldovich dynamics. 
Therefore, the conditions of convergence of various integrals 
are the same as in Sect. 17.11 

We now consider the matter bispectrum. At tree-order, 
the diagram (a) of Fig. [H] gives back the diagram (a) of 
Fig. El Thus, the explicit computation of Ea. (p5)) recovers 
Ea. ((9T|) for the Zeldovich dynamics. At one-loop order, we 
no longer have such a one-to-one correspondence between 
both expansions since we have eight one-loop diagrams in 
Fig. [6] and only four one- loop diagrams in Fig. [3] First, let 
us consider the contributions which have the form of a triple 
integration over the linear power spectrum, as in Ea. (j95p 
associated with diagram (e) of Fig. [S] One can see at once 
that only diagram (f) of Fig. [H] has this form and its explicit 
calculation gives back Eg. ([55)1 . Next, contributions with a 
double integration, as in Eg. ([94)) . arise from diagrams (b) 
and (g) which read as 

S(^) - -PL(fci) ydqidq2 Soim + q2 - ka) x 
2(ki.k3)(k2.qi)(k2.q2) + (2fc2(k2.k3)-fci(ki.k2))(qi.q2) 



PL{qi)PL{q2) + 5 perm. 



(k2.qi)(k2.q2) 

91 '72 



(100) 



i(fc) = 



dlnP^o(fc) 
dlnfc ' 



(97) 



we see that the one-loop contributions to the power spec- 
trum and the bispectrum converge if n(0) > —1 and 
n{oo) < —1. However, as is well-known (Vishniac 1983; 
Jain & Bertschinger 1996), the infrared divergences at 
q —^ cancel out and the sum of all contributions is fi- 
nite for n(0) > —3. This can be easily checked from the 
explicit expressions (|90 p - ([M)) . On the other hand, the time- 
dependence of various contributions is 



P"™ oc Pi^cr, 



B 



8' 



tree ^ n4„4 
OC 1^ (Tg , 



plloop ^ ^4^4^ 



plloop ^ ^6 fi 



(98) 
(99) 



where erg is the normalization of the linear power spectrum 
today at 2 = (defined as usual from the rms linear density 
contrast within a sphere of radius 8/i~^Mpc). Thus, we have 
the scaling P^ oc D^a^ for a fixed shape of the linear power 
spectrum. Higher-order terms scale as higher powers of P^ 
hence of D^ct^. 



and 



P^^') --PL(fci) y"dqidq2 Soim + q2 - ^2)PL{qi)PL{q2) 
"2fc2g2(k3.q2) + 2(ki.qi)(k3.q2)2 



+ 



6kfqfq1 
(ki.qi)[-fc|g2 + (fc2^ 



2fc| 



ql 



<?l)(k3.q2) 



Gkfqfqi 



(k2.qi)(k2.q2) , _ 
X ^-^5 h 5 perm. 

qiq2 



(101) 



Then, summing P^^' -f B^^^ and writing the result in a 
symmetric form over {qi,q2} we recover the contribution 
of Eg. ([M| . Finally, contributions with a single integration, 
as in Eg. fM]) . arise from diagrams (c), (d), (e), (h) and (i), 
which read as 



S(c) 



D ii. \u ii. s (ki:k2)(kijc3) / , p / N 
-PL{k2)PL{k3) T7^72 /dqPL(g) 



(k2.q)2 + (k3.q) 



2 



2 perm., 



(102) 
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3fc| fc| + 2 (fc2 + kl ) (k2 .kg ) + (k2 -ka )2 and 



X JdciPLiq) 



1 2 



(k2.q)2 + (kg.q) 



2 perm. 



(103) 



B^'^^ = -PL{k2)PL{k3) 



2 A^2 ^3 



^i)(k2.k3) 



dqPL(g) 



(ki.q)2 



BW=Pi(fc2)Pi(fc3) JdqPLiq) 



2 perm., 



(k2.q)(k3.q) 

1 2 ^2 



(104) 



,(kl.q)^ 



kfq^ + Sk^ki + 2{ki + fc|)(k2.k3) + (k2.k3) 



21 



3(72(ki.q)[fc|(k2.q) + fc|(k3.q)] 



+2 perm., 



(105) 



and 



"(ki.q)(k2.q)fci(2fci + 3g2) 



BW^P^(fc2)Pi(fc3) y'dqP^(g)[ ^^^^ 

, (ki .q) (k3 .q)fci {2kl + 3g2) - (ki .q)^ (k2 .k3) (fc2+ ^32) 



12/^2 k^ 

(k2.q)(k3.q)[2fc2g2_3(k,.q)2]- 



2 perm., 
— q and k2 



(106) 
ks, and 



1 2 k'2 k^ Q 

where we used the symmetries q 
the property ki -I- k2 -I- k3 = 0, to simplify the expressions. 
Then, we can check again that the sum 
pCi) of Eqs. pfe)) - pU51) gives back Eq.® obtained 

from the standard perturbative expansion. This consistency 
check allows us to check these expressions obtained through 
two different methods and to track possible mistakes in the 
computations. 

We again find that all one-loop contributions to the bis- 
pectrum converge if n{0) > —1 and 71(00) < —1, except for 
B*^''-' and i?^*-' which require n{oo) < —3. As for the stan- 
dard perturbative expansion (which we recover by summing 
all terms) infrared divergences compensate so that we only 
need n(0) > —3. On the other hand, the UV divergences of 
B'^^^ and B*-'^ also compensate so that we recover the stan- 
dard constraint n(oo) < —1. Thus, although the standard 
expansion of Sect. 17.11 and the expansion over of this 
section are actually identical, the contributions obtained 
at each order are split in different ways between different 
diagrams so that the convergence properties of individual 
parts can be different for the two schemes. In particular, 
we have seen that for the Zeldovich dynamics one obtains 
an artificial UV divergence for —3 < n(oo) < —1 which dis- 
appears by summing all diagrams. For numerical purposes, 
since the contributions (h) and (i) of Eas. (|105p - (|106[) are 
not well-defined because of this UV divergence, we intro- 



duce the regularized contributions B^^^ and B''^' . They are 
equal to B^'*) and i?'-*' from which we subtract the UV di- 
vergent part, that is, we remove the part that scales as q'^ 
from the terms in the brackets in Eas. (|105p - (|106p . This only 
leaves terms that scale as q~^ and gives 



B 



(h) 



'PL{k2)PL{k3 



Zklkl + 2(A:|-Hfc2)(k2.k3) + (k2.k3 



6/^2 k^ 



IdqPUq) 



(k2.q)(k3.q) 



2 perm.. 



(107) 



B^j^'^ = -PL{k2)PL{k3) 

X JdqPLiq)^^^^ + 2 perm 



2fc|fc| + (k2.k3)(fc2 + fc2 + fc2) 



1 '2iki2^ /c^ 



(108) 



Since the UV divergences of S'^'*' and i?^*) compensate we 



have pC') + = B^^' + B^^ and we can use B^^' and 
B^^^ instead of B^') and p('). 



7.3. Direct steepest-descent method 

We now consider the direct steepest-descent method pre- 
sented in Sect. 15.11 First, the two-point functions C and 
R are obtained from the diagrams of Figs. [5][2] that corre- 
spond to Eas. ([7T|) - ((73)) . Thus, in terms of nonlinear two- 
point functions we have only one one-loop diagram for ei- 
ther C or i?, which corresponds to an infinite series over the 
linear two-point functions Cl and Rl- The results obtained 
at one-loop order have been described in detail in Valageas 
(2007b). Thus, the nonlinear response R reads as 



R{xi,X2) = Rl{xi,X2) cos[uj{k){Di - D2)], 



(109) 



where Rl is the linear response given in Eg. p^ . D = 
is the linear growth factor and uj(k) (where fc = fci ~ ^2 is 
the wavenumber) is given by 



uj(k) — k(7y. 



(110) 



Here a'^ is the variance of the one-dimensional linear dis- 
placement field slo given in Eg. ([55)1 . Note that the nonlin- 
ear response R only depends on the linear power spectrum 
through this velocity variance. From the expression of the 
linear response Rl, Ea. (|109p gives 



R{k;Di,D2) = e{Di-D2)cos[ij{k){Di-D2) 



2 

^ V^'2 

£=1 ^ 



R: 



with 



ai 



1, a2 = 0, 



(111) 



(112) 



and the two constant matrices and ^(2)5 associated 
with the linear growing and decaying (here constant) modes 
are 



^(1) = ( 1 



and i?(2) = 



1 -1 




(113) 



The expression of the two-point correlation C is more 
intricate and it depends on the details of the linear power 
spectrum. From Eq. (|109p the first term of Eg. ([75)1 reads 
exactly as (Valageas 2007b) 

R X Cl{vi) X R^ ^ Cl{xi,X2) cos(cjDi) cos{ujD2). (114) 

The second term of Eq. ([75)1 has a more intricate structure 
but we shall use the qualitative approximation 



R.n.R^ ~ CL(xi,X2)sin(wi:>i)sin(cji:>2). 



(115) 



As seen in Valageas (2007b), the exact one- loop result con- 
tains an additional term but it has no qualitative effect on 
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C hence we shall use the simple approximation obtained power spectrum and on the wavenumbers kj through the 



from the sum of Eqs. (lll4p - (lll5p . 

C(xi,X2) ~ CL{xi,X2)cos[uj{k){Di 
which gives: 



three combinations kjU^. The explicit expressions of T, 



(a) 



(116) 



are 

(a)_ 



{I) 



T, 



(1) 



sinc[_D(— oji + UJ2 + ^3)] + sinc[D(a;i — UJ2 + 1^3) 



C(k;Di,D-,)^DiD2Cos\Lu(k)(Di-D2)]PLo{k)( ] ]].(n7) , ■ 1 m 1 • rn/ , , m\ no1^ 

^ ' ^ L V A i ^/j ^yi ij ' +smc[L'(tJi + - wsjj + smc[L'(wi -I- W2 + wajj ^ (121) 



Equations (I109p . (lll6p . also correspond to a simple effec- 
tive dynamics where the nonlinearity is reduced to a simple 
non-Gaussian random advection by a large-scale flow which 
gives rise to the cosine factors (Valageas 2007a, b). Then, 
these cosine factors simply describe the (weak) decorrela- 
tion, or loss of memory, due to this random advection (com- 
pare with Eqs. ([5^ . ([55)1 V Note that the time dependence 
obtained from the steepest-descent method is quite differ- 
ent from the one obtained in Eq. ((98l) since Eas. (lll4p - (|115p 
give for the equal-time power spectrum 



(118) 



where we label as the "tree" contribution the first diagram 
of Fig. [21 associated with Eq. (|114p . and as the "Hoop" 
contribution the second diagram of Fig. [51 associated with 
Eq. (|115p . Therefore, we now find that both contributions 
exhibit the same scaling in the nonlinear regime and they 
only grow as -D^u| instead of D'^a'^ as in Eq. ([55)) . 

The three-point correlation is now given up to one- 
loop order by the five diagrams of Fig. [TUl As explained 
in Sect. 15.11 this simplification with respect to the expan- 
sion over powers of Ks, studied in Sects. l4l and 17. 2[ comes 
from the fact that one-loop diagrams associated with a 
"renormalization" of the two-point functions C and R are 
included in the tree-diagram (a) of Fig. [10] and only the 
four diagrams associated with the "renormalization" of the 
three-point vertex Kg are new structures with respect to 
the previous order (tree-order). Let us first consider the 
tree-diagram (a). Its explicit expression is of the form of 
Eq. (|65p but using the nonlinear two-point functions R and 
C. From Eqs. (|llip . (|117p . we see that the integral over time 
rj'i in Eq. ([65|l gives rise to the two quantities (for £ = 1,2): 



(119) 



This defines the two time-dependent functions T^^^^^ of 

the three variables LOjD, associated with diagram (a). In 
Eq. (|119p we considered the case of equal-time statistics 
Di ^ D2 = D3 = D and we introduced: 



UJj = Lj{kj) = kjUy. 



We also normalized T^^} so that for 



(120) 

uiiuu i^^i = we have 

0) = 1. Thus, by putting T^^^ = 1 we must re- 
cover the result ([OTjl obtained by integrating Eq. ([S5)) over 
the linear two-point functions, as in the expansion scheme 
of Sect. 17.21 where we expanded over powers of K^- This 
provides a convenient consistency check. Note that all the 
time-dependence is included in the standard prefactor 
and the function T^^^^ , which only depends on the linear 



and 



T^2) = ^|C0SC[L'(-LJ1 +UJ2+ t^3)] + COSc[D(cJi - UJ2+ ^3)] 

+cosc[D{uJi + Uii2— iJ-'3)] + cosc[Z)(a;i + UJ2 + ^3)] > , (122) 



where we introduced the functions 

, , sin(x) , , 2(1 — cos(a;)) . 9, a;, 

sinc(x) = — cosc(a; =-i = sinc^ ( - ) . (123 

X a;^ 2 

We can check that t'Ij^ (0) = 1 , and in the highly nonlinear 
regime we have 



(124) 



(with I ^ ±Wj| ^ Lo). Then, diagram (a) of Fig. [10] reads as 



B(-)=Pi(fc2)Fi(fc3) 
-\-2 perm. 



Tf)fc?(k2.k3) + rg[A;|fc|-(k2.k3)2] 



"-2 '^3 



(125) 



We can check that by substituting t'"^^'^^ = 1 (and using ki -|- 

k2+k3 ~ 0) we recover the standard tree-order contribution 
(I9ip . Moreover, in the nonlinear regime we have the scaling 



(126) 



since u oc erg. Thus, the nonlinear corrections resummed 
in the large-A'' steepest-descent method have modified the 
scaling over D of the tree diagram (a) as compared with 
the standard perturbation theory result ([9T|) which scales 
as the first term in ([M)) . 

Next, let us consider the first one-loop diagram (f) of 
Fig. [TO] We can again factorize the time-dependent parts 
which read as 



^ ' dD' ( dV' 



^ 



U 



~ (2-a.J(2-a,J(2-a.3)^(^-^-^3)('-.^)- (^27) 
i f) 

This defines the functions of the three variables cojD. 

They are parameterized by the indices £j of the response 
matrices R(ej) of Eq. (|113p . Since the diagram (f) involves 
three response functions, there are 2^^ = 8 such functions 
T^^Py However, thanks to the symmetry of diagram (f), T^^P^ 
can be split into a product of three terms as 

TtP,i,,i,)i^,D)=TlP^{uj,D)TlP^^^^^ (128) 
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with 



tI^^\^D) = (2 - a^) ^x^'"' cos[ujD{l - x)]. (129) 
The integration over x yields 

tI^(^{u;D) = sinc(wi:>), T^J^\ljD) = cosc(wD), (130) 
which gives the nonlinear scaling 



(131) 



and higher powers of 1/{Dlu) arise if some indices £j are 
equal to 2. Then, diagram (f) reads as 



= ydqidq2dq3 Soiqi + q2 - ki)^z3(qi -)- qs + k2 



Tll^\LU,D)kf{q,.q,)^ 




(qi.q2)'] 




9 9 

mi 






-T[l^{u2D)[qlql- 


(qi.qa)'] 




9 9 

mi 




-T(//(u;3i?)fci(q2-q3) 


+ Tll\u,D)[qhl 


-(qi.qa)'] 


mi 



(132) 



We can check that by substituting T^'^^^ = 1 we again 

recover the corresponding one-loop diagram (f) of the 
expansion over Kg of Sect. \7.2\ which is also given by 
Ea. (|95|) obtained from the standard perturbative expan- 
sion of Sect. 17.11 Then, we obtain for this one- loop diagram 
(f) the scaling 



DiUj > 1 : B'-f^ 



D'4 



3^3 



(133) 



Thus, the scaling of this one-loop diagram is again modi- 
fied as compared with the corresponding one-loop diagram 
obtained in the standard perturbation theory, which was 
equal to Eg. ([55)1 and which scaled as the second term of 
(|99|) . Moreover, contrary to the standard expansion scheme, 
the one-loop diagram shows the same scaling as the tree di- 
agram (compare with Eq. (|126p ). 

The three other one-loop diagrams of Fig. [10] can be 
computed in the same fashion. In particular, because of 
the triple integration over inner times D'p weighted by the 
factors cos[wj(D — D'^)] as in Eq. (|127p . we can check that 
all time-dependent factors T obey the same scaling as in 
Eq.dni]), 



Dujj > 1 : 



1 



(134) 



with possible higher powers if some Ij are equal to 2. 
Therefore, all one- loop diagrams scale as Ea. (|133p and as 
the tree diagram of Ea. (|126p computed at this same one- 
loop order. Thus, we recover the nice improvement over the 
standard perturbation theory already obtained for the two- 
point correlation: the one-loop correction does not grow 
with an additional power as compared with the tree- 
order result and seems better behaved. Unfortunately, as 
shown in Valageas (2007b), at higher orders exponentially 



growing terms appear within this steepest-descent expan- 
sion scheme. Therefore, at two-loop order and beyond, this 
large-iV method becomes even more ill-behaved than the 
standard expansion in the sense of exhibiting strongly grow- 
ing high-order terms. 

On the other hand, we can note that the growth as I?^cr| 
of the two-point correlation (see Eq. (|118[) ) and as -D^cr| of 
the three-point correlation breaks physical constraints, that 
is, it cannot be achieved by any matter density distribution. 
Indeed, the positivity of the matter density p(x) implies 
that C53 should grow at least as (C52)^ in the highly non- 
linear regime, see Eq.(6) in Valageas (1999). This implies 
that the strong constraint p((x) > is not satisfied by the 
large- steepest-descent expansion (nor by the usual per- 
turbative expansion). The growth factors and are 
rather reminiscent of generic unconstrained fields '0 with 
a time-dependence which would approximately factorize as 

Finally, we must take care of the UV divergences of di- 
agrams (h) and (i) already encountered in Sect. 17.21 In the 
expansion over Ks, described in Sect. 17.21 the UV diver- 
gences of diagrams (h) and (i) exactly compensate and we 

could as well use the regularized contributions B^j^^ and 

b'^j^' of Eqs. pHTIl - piIg)) . For the steepest-descent scheme 
this is no longer true because the divergent parts are mul- 
tiplied by the independent functions T*^^) and T^'^. On the 
other hand, at order Pf^ we could substitute T^^^ = T'*) — 1 
into diagrams (h) and (i), which already involve a product 
over three terms P^, so that the divergences again com- 
pensate. In fact, at order P| we must recover the results of 
the standard expansion scheme hence we know that such a 
UV divergence can only appear at higher orders (if at all) . 
Therefore, in a fashion similar to the procedure of Sect. 17. 2[ 

we define regularized contributions -B^^ and -B^'' by sub- 
tracting from B^'*) and S*-'^ their UV divergent part. This 
makes no difference for the sum of contributions (h) and 
(i) at order (but we have removed higher-order terms 
which happened to diverge and did not cancel in the sum). 
This gives: 



B 



_ Tl'^l^2{\.Mkl + Tl^l^Smki - (k^.ka)^] 



'PLik2)PLik3) JdqPL{q) 



(k2.q)(k3.q) 



2 perm. (135) 



and 



T«,)(k2.k3)fc2+T, 



(h) 
(212) 



(k2.k3 



6 kn kri 



xPL(fc2)PL(fc3) dqPLiq) 



1,2 n,3 

,(kl.q) = 



2 perm. 



(136) 



Of course, we can check that by substituting T^') = T^') = 
1 into Eqs. (|135|) - (|136p . and using ki + k2 -f k3 = 0, we re- 
cover the regularized contributions (I107p - (jl08p introduced 
in Sect. 17.21 for the expansion over Kg. On the other hand, 
this analysis shows that the partial resummations involved 
in the large- iV methods have introduced spurious UV di- 
vergences. For the Zeldovich dynamics, we know that such 
UV divergences are not seen in the exact nonlinear results 
that can be obtained from the solution l|35p of the dynam- 
ics. Therefore, they must cancel out by including further 
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diagrams or by expanding back over powers of Pl, as in 
the standard perturbation theory. The fact that they do 
not cancel out automatically may be seen as a shortcoming 
of resummations schemes such as these large- TV methods. 
Indeed, the use of partial resummations is intended to im- 
prove the behavior of the predicted nonlinear correlations as 
compared with the results of the standard perturbative ex- 
pansions. For instance, one would like to recover the damp- 
ing of the response function in the nonlinear regime and a 
better control of higher-order terms. At one-loop order this 
is partly achieved for the bispectrum as all contributions 
show the same scaling (|133p instead of the hierarchy ((99)) . 
However, this comes at the price of an increased sensitivity 
to the smallest scales that even leads to a divergence for 
n{oo) > —3. Nevertheless, since this divergence disappears 
in the expansion over it is always possible to remove it 
by subtraction as in Eas. (|135p - (|136p . 



8. Gravitational dynamics 



As seen in Sects. [2][5l we can apply to the gravitational 
dynamics the expansion schemes described in Sect.[7]for the 
Zeldovich dynamics. We do not give the explicit expressions 
of various diagrams here as they are rather long, except for 
the expansion over Kg in the Appendix. 

8.1. Standard perturbation theory 

For the standard perturbation theory recalled in Sect. [31 the 
kernels Fn are computed from recursion relations derived 
from the equation of motion ([8]). Following the procedure 
described in Sect. [31 this gives the well-known results for 
the matter power spectrum and bispectrum that we do not 
repeat here (Bernardeau et al. 2002 and references therein; 
Scoccimarro & Frieman 1996a, b; Scoccimarro 1997). The 
conditions of convergence of individual diagrams are the 
same as for the Zeldovich dynamics of Sect. 17. li that is, 
n(0) > — 1 and n(oo) < — 1. Again, since the Galilean in- 
variance applies to both dynamics, the infrared divergences 
cancel out as we sum all diagrams so that we only re- 
quire n(0) > —3 (Vishniac 1983; Jain & Bertschinger 1996). 
Finally, the dependence on the linear growth factor D{ri) 
and on the amplitude as of the linear power spectrum are 
as in Eqs. ([551) - ([M1) . We do not plot the results here since 
they are identical to the results obtained from the expan- 
sion over Kg displayed in the next section. 

8.2. Expansion over cubic interaction 

The expansion over powers of Ks proceeds exactly as for 
the Zeldovich dynamics studied in Sect. 17.2] For the power 
spectrum there is again a one-to-one correspondence be- 
tween the diagrams (b) and (c) of Fig. \5\ and Fig. ^ For 
the bispectrum, the correspondence is again between the 
groups of diagrams {{a),{b + c),(d),(e)} of Fig. [3] and 
{(a), ic + d + e + h + i),{b + g), (/)} of Fig. E For com- 
pleteness, we give the expressions of these diagrams in the 
Appendix. 

These results are similar to those obtained for the 
Zeldovich dynamics in Sect. 17.21 except that the ratio- 
nal functions that appear in the integrals are more intri- 
cate. In particular, the denominators do not factorize as 
products of kj and qj so that angular integrations require 




k [h Mpc-i] 

Fig. 15. The reduced bispectrum Q defined at one-loop or- 
der by the rational function (jl39p . We show Q as a function 
of wavenumber k for equilateral triangles fci = ^2 = fcs = k. 
The data points are taken from Fosalba et al. (2005), for 
a ACDM simulation at z = 0. The letters correspond to 
the various contributions, such as B^"' /{S^'"" + 5'"°°?), as- 
sociated with tree and one-loop diagrams obtained for the 
bispectrum B within the expansion over Kg . The solid line 
Q is the full reduced bispectrum Q (i.e. the sum of all these 
contributions). It is also equal to the prediction of the stan- 
dard perturbation theory. The rising dotted line Qpcs is the 
phenomenological model of Pan et al. (2007) . 



more care. The conditions of convergence are the same 
as for the Zeldovich dynamics, that is, n(0) > —1 and 
n(oo) < —1, except for contributions i?^''-' and B*^*-* which 
require n(oo) < —3. Again, the infrared and ultraviolet di- 
vergences compensate so that the sum is well-defined for 
n(0) > -3 and n{oo) < -1. As in Eqs. ([TU7|) - piI5)) . it is 

convenient to introduce the regularized contributions sj^'^ 

and sj^' by subtracting from S'''^ and S'*' their UV di- 
vergent part. 

From the bispectrum B{ki, k2, k^) and the power spec- 
trum P(fc) we define as usual the reduced bispectrum 
Q{ki, k2, fcs) by 



Q{ki, fc2, fcs 



B(fci,fc2,fc3 



^'(^2)^(^3) + 2 perm. 



(137) 



From Eq. (|5ip . in the quasi- linear regime it is independent of 
time and of the normalization of the linear power spectrum. 
Let us note S{ki, k2, k^) the denominator of Eq. (|137p : 



5(fci, fc2, ks) = P{k2)P{k3) + 2 perm. 



(138) 



Within the standard perturbative expansion over powers 
of Pl , or the equivalent expansion over Kg , the ratio (|137p 
reads up to one-loop order as 



Q 



^trcc _j_ ^lloop 
gtvcc ^ 5'lloop 

^ticc ^lloop ^trcc ^lloop 



^tioc ^tr 



(139) 
(140) 
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0. 1 1 

k [h Mpc-i] 

Fig. 16. The reduced bispectrum Q at one-loop order for 
equilateral triangles as in Fig. 1151 but computed from the 
expansion p40p instead of the ratio (|139p . The symbols are 
as in Fig. [15] but the contributions to the bispectrum B as- 
sociated with tree and one-loop diagrams are now divided 
by S^'^'^^, as in the first two terms of Ea. (|140p . instead of 
^^tree ^ ^iioop)^ g^^^^ ^^le ucw curvc labeled -BS/S^ corre- 
sponds to the last term of Ea. (|140p . 



where in the second line we have expanded the ratio B/S of 
the first line. The various contributions to the bispectrum 
have been described in the previous sections, whereas the 
denominator S is given up to one-loop order by 

^t'- = Pi(fc2)Pi(fc3) + 2perm. (141) 
^iioop ^ p^(A;2)P"°°P(A:3) + 5perm. (142) 

We show in Fig. [12] the results obtained for the reduced 
bispectrum Q as defined by Eq. (|139p . that is, we keep the 
ratio (|139p without expanding again as in Eq. (|140p . We plot 
Q as a function of wavenumber k for equilateral triangles 
ki = k2 = = k. The data points are taken from Fosalba 
et al. (2005). They correspond to a ACDM simulation from 
the Virgo archive at redshift z The four dashed curves 
labelled {a,c + d + e + h + i,h + g, f} are the contributions 
of the various diagrams obtained for B in the expansion 
over Kg- Diagrams of the same form (i.e., that involve the 
same number of integrations over the linear power spec- 
trum) have been gathered and divided by the denominator 
^tree _|_ ^Uoop^ Their sum is the full reduced bispectrum 
Q shown by the solid line. The rising dotted line Qpcs 
is the phenomenological model of Pan et al. (2007). It is 
based on the scale transformation introduced by Hamilton 
et al. (1991) and it uses as input the nonlinear power spec- 
trum, obtained from the fit to numerical simulations given 
by Smith et al.(2003). We show in Fig. [16] the result s ob- 
tained for the same case when we use the expansion (|140p 
instead of the ratio (|139p . The new curve labeled —BS/S^ 
is the last term of Eq. (|140[) . 

Let us recall here that the predictions for the bispectrum 
B, hence for the reduced bispectrum Q, obtained from this 
expansion over Ks , are equal to the results of the standard 
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Fig. 17. The reduced bispectrum Q as a function of the 
angle 6 between wavenumbers ki and k2. Panels (a), (b), 
and (d) have k2 = 2ki whereas panel (c) has ^2 = ki. We 
probe deeper into the nonlinear regime by going from panel 
(a) to panel (d). The solid line Qi corresponds to the ratio 
(|139p whereas the solid line Q2 corresponds to its expansion 
(|140p . They are almost identical in panel (a). The dashed 
line is the tree-order result (|143p . The dotted line Qpcs 
is the phenomenological model of Pan et al. (2007). Data 
points are taken from Scoccimarro & Couchman (2001), 
for panels (a) and (b), and from Fosalba et al. (2005), for 
panels (c) and (d). They correspond to two different ACDM 
cosmologies. 



perturbation theory. Note that the tree-order result, given 
by 

otree 



also corresponds to the horizontal line (a) in Fig. [TH] It 
is constant since both iJ^rco 5'*'''=<= scale as PL{k)^. We 
can check that the reduced bispectra Q obtained from the 
two Eqs. (|139p - (|140p agree on quasi-linear scales, up to 0.2h 
Mpc-\ where A|(fc) ~ 1 (here A|(fc) = A-nk^P{k) is the 
power per logarithmic wavenumber). At smaller scales the 
ratio (|139p goes to a constant (that is, a value that is inde- 
pendent of the normalization of Pl but shows a weak de- 
pendence on the shape of the linear power spectrum, since 
both g" °°P a nd 5"°°p scale as P|), whereas the expanded 
version (|140p grows as A^. This leads to a better agree- 
ment of Eq. (|140p with the numerical results up to 0.4/i 
Mpc^^, where A|(fc) ~ 4. However, it is clear that this 
better match has no strong foundations. We can also note 
that in the nonlinear regime there are some cancellations 
between various contributions, for both Eqs. (|139p - p40p . 
On the other hand, it appears that the phenomenological 
model of Pan et al. (2007) works best, up to 0.8ft. Mpc~^, 
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where A|(fc) ~ 17. However, at fc ~ 0.1 - 0.2h Mpc^^ it 
may slightly overestimate Q. 

We show in Fig.[T7]the reduced bispectrum Q as a func- 
tion of the angle 9 between the two wavenumbers ki and 
k2: 



cos(6') = 



ki.k2 

fclfc2 



fcg = Jkl + kl + 2fcifc2 cos(6'). (144) 



In the lower left panel (c) we have fc2 — ki whereas in 
other panels ^2 — 2ki . We probe smaller scales farther into 
the nonlinear regime going from panel (a) to panel (d). 
On large scales, all curves match the tree-order result and 
numerical data, as seen in panel (a), and show a strong de- 
pendence on the angle 9. On smaller scales, the dependence 
on 9 decreases, as shown by the simulations. This behavior 
is captured by the phenomenological model of Pan et al. 
(2007) but not by the one-loop predictions of Eqs. (|139p - 
(|140p . In agreement with Figs. [T51 - [TBI the ratio (I139P yields 
a smaller value for Q than the expansion (|140p . but none 
of these two predictions matches the data for A|^(fc) > 0.3. 



The results displayed in Figs. 



agree with previous 



studies of the standard perturbation theory, discussed for 
instance in great details in the review of Bernardeau et al. 
(2002). 



8.3. Direct steepest-descent method 

For the gravitational dynamics, the nonlinear response 
function is no longer given by Eq. (jl09p at one-loop order 
in the steepest-descent expansion. However, Eas. (|109p and 
(|117p still provide a reasonable approximation of the exact 
one-loop result, as seen in Valageas (2007a). Thus, in this 
article we use the approximations (|109p and (|117p to in- 
vestigate the behavior of the steepest-descent method for 
the gravitational dynamics as well. Then, from Eq. fTO)) the 
powers and the matrices of Eqs. (|112p - (|113p must 
be replaced by 



ai 



a2 



3 2 
3 2 



R 



(2) 



2 

-3 



(145) 



(146) 



Next, the computation of various diagrams proceeds as in 
Sect. 17.31 In particular, since the linear growing modes are 
identical for both dynamics (exponent ai — 1), the time- 
dependent functions T^j^^^ are identical if all £j are equal 

to 1. Then, we recover the scalings (|118p for the power 
spectrum and (|126p . (|133p . for the bispectrum. 

Again, in order to handle the UV divergences associated 
with diagrams (h) and (i) it is necessary to introduce the 

regularized contributions b\^^ and , as in Eqs. (|135p - 
(|136p . As discussed in Sect. 17.31 this makes no difference at 
order as the UV divergence is associated with higher or- 
der terms. On the other hand, for the gravitational dynam- 
ics this UV divergence may have a physical meaning. That 
is, it could be a signature of the breakdown of the equations 
of motion beyond shell-crossing and may not compensate 
at higher orders. Indeed, as noticed in Valageas (2002), 
within the standard perturbative expansion one also en- 
counters UV-divergent diagrams beyond a finite order if 
n(oo) > —3. Since the large- expansion has performed a 



partial resummation of an infinite series of diagrams it is 
not too surprising to encounter such UV-divergent terms. 
For the Zeldovich dynamics this is not the case, since we 
know from the exact nonlinear results that these diver- 
gences must eventually compensate once we include all 
diagrams of a given order, even though the Eulerian de- 
scription also breaks down beyond shell-crossing. However, 
for the gravitational dynamics, where there is no such ex- 
plicit nonlinear solution, it is not known whether these 
divergences must also fully cancel out (nevertheless they 
must compensate at least up to the last regular order pre- 
dicted by the standard expansion scheme) . We do not give 
the explicit expressions of the diagrams obtained in this 
fashion here as they are too lengthy. They are similar to 
those obtained within the expansion over Ks, given in the 
Appendix, with the addition of time-dependent functions 
T^^'^'^'^'^'^\ijjjD), as described in details in Sect. 17.31 for the 
simpler case of the Zeldovich dynamics. 

Since we use Eqs. (|109l) and pi7p for the nonlinear 
two-point functions the reduced bispectrum Q defined in 
Eq. p37|) reads as 



(5cos(fcl,fe,^3) = 



gcos(fcl,fe2,fc3) 

PL(fc2)/'L(fc3) + 2perm.' 



(147) 



where the subscript "cos" stands for the steepest-descent 
method with the original cosine cutoffs of Eqs. (|109p . (|117p . 
At this one-loop order, Bcos is given by the diagrams of 
Fig. [TUl as discussed above. 



C3^ 



- 




0.1 



k [h Mpc-1] 



Fig. 18. The reduced bispectrum Q function of 

wavcnumber k for equilateral triangles ki — k2 = k^ — k, as 
in Fig. [ini We show the prediction of the one-loop steepest- 
descent method from Eq. p47p (solid line Qcos) and the re- 
sults obtained by using a Gaussian cutoff as in Eqs. ((82|) . 
(dotted line Qcauss)- 



As noticed in Sect. [U it is possible to consider an al- 
ternative resummation scheme for the bispectrum, based 
on the diagrams of Fig. [TU] associated with the large- A^ 
steepest-descent method, but using the two-point functions 
((82| . ((86|) . given by the large-fc partial resummation asso- 
ciated with the random uniform advection of large-scale 
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Fig. 19. The reduced bispectrum Q as a function of the 
angle 9 between wavenumbers ki and k2, as in Fig. [TOl 
The curves Qcos (solid line) and Qcauss (dotted line) are 
the predictions obtained from the one-loop steepest-descent 
method by using, i) the original nonlinear two-point func- 
tions with a cosine cutoff, ii) the two-point functions with 
a Gaussian cutoff. 



structures. This amounts to replace the cosine cutoffs of 
the one-loop steepest-descent predictions (|109p . (|117p . by 
the Gaussian cutoffs of Eqs. ([5^ . ([M)) . Therefore, the struc- 
ture of the various contributions to the bispectrum B is the 
same as for the genuine steepest-descent method discussed 
above. The only change is that the time-dependent func- 

are different as they involve integrals 

over Gaussian factors instead of cosines. Nevertheless, the 
scalings (jl24l) and (|134p remain valid. 



DiUi > 1 : T 



(tree) 



1 



^ rTi(lloop) 



where we noted by t(*''°°^ the time-dependent functions as- 
sociated with the tree diagram of Fig. \W\ and by t("°°p) 
the time-dependent functions associated with the one-loop 
diagrams of Fig.fTUl Indeed, the scalings (|148p simply come 
from the single and triple integrations over past times asso- 
ciated with the tree and one-loop diagrams. For large Doj, 
only recent times D' , with {D — D') <i l/w, contribute 
because of the Gaussian or cosine cutoffs. This leads to 
Eq. (|148p . Therefore, we again recover the same scalings 
(|126p and (|133p for the tree and one-loop contributions to 
the bispectrum B. Then, the reduced bispectrum Q is again 
given by Eq. (|147p . but using a subscript "Gauss" to denote 
the Gaussian cutoffs associated with the resummation of 
Sect, m 

Thus, for the equal-time bispectrum the form of the 
different-time decay of two-point functions does not really 
matter since it merely yields the scalings (|148p . Of course. 



for the different-time three-point correlation we recover the 
decay of the two-point functions. For instance, we obtain 
up to one-loop order from the diagrams of Fig. [10] 

Djtoj > 1 : B ~ e-5[(^2--Di)'-'2'+(^3-i3i)^^il + 2perm. (149) 

if we use the Gaussian decay. 

We show in Fig. [18] the results obtained from the two 
methods described above for the reduced bispectrum as a 
function of wavenumber k for equilateral triangles. We can 
see that the the curves agree on large scales with each other 
and with the numerical simulations, up to 0.2ft, Mpc~^. The 
comparison with Figs. I15[ 1161 shows that the approxima- 
tions (fT09)) . (fTTT)) . or (ISll), (US]), for the nonhnear two-point 
functions, are accurate enough to recover the deviations 
from the constant tree-order result associated with the one- 
loop order corrections. On small scales, the reduced bispec- 
tra Qcos and Qcauss show a slow decay, in agreement with 
the scahngs ([TT8| and (fT26ll . (fT33| . which give 



Dujj > 1 : Q 



1 



Da. 



(150) 



The curve Qcos displays some oscillations due to the co- 
sine factors whereas the curve Q Gauss shows a monotonous 
decay due to the Gaussian factors. 

We show in Fig. [TO] our results for the bispectrum as a 
function of the angle 6 between wavenumbers ki , k2 . We 
can see that the results are similar to those obtained in 
Fig. [T7] from the expansion over Kg (equivalent to the stan- 
dard perturbation theory). In particular, they do not man- 
age to reproduce the flattening of the reduced bispectrum 
as a function of 9 on smaller scales. Therefore, the pre- 
dictions obtained in this manner from the steepest-descent 
formalism do not bring a significant improvement over the 
standard perturbation theory, as shown by the comparison 
with Figs. [IMH] 

8.4. Using the nonlinear power spectrum 

We do not compute numerically in this paper the predic- 
tions associated with the 2PI effective action method of 
Sect. 15.21 because there is no simple approximation such 
as Eqs.([l09l), (fTTTp . or Eqs.([82l), for the two-point 

functions. Moreover, the integrals over time no longer fac- 



(148) torize through functions T, 



{ujjD), since the inter- 



nal lines of the one-loop diagrams are nonlinear two-point 
functions. Thus, even if the latter are approximated by 
Eqs.HnU), (fTTT)) . or Eqs.((82l), one encounters integrals 
over internal times D'^ and wavenumbers qj of terms such as 
Q-iDj-D'-)q^cr^,/2 instead of power-law prefactors (D'^/Dj)". 
All this would make the computation of the 2PI effective 
action predictions much more complex. 

However, to investigate the effect of using more accu- 
rate nonlinear two-point functions in all lines of the di- 
agrams of Fig. 1101 we modify the approximation Qcauss 
discussed above by substituting the fit Pnl of the nonlin- 
ear power spectrum, obtained from numerical simulations 
by Smith et al. (2003), into the diagrams. In other words, 
we use the exact equal-time nonlinear two-point functions 
in all lines, and the dependence on the time-difference on 
external lines is given by the Gaussian factor of Eqs. ([5^ . 
(|86)) . We keep the power-law time factors in the internal 
lines so as to keep the factorization through the functions 
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0. 1 



k [h Mpc-i] 



Fig. 20. The reduced bispectrum Q as a function of 
wavenumber k for equilateral triangles ki = k2 = = k, 
as in Fig. 1181 but using the nonlinear power spectrum from 
Smith et al. (2003) for the two-point correlations within 
the one-loop steepest-descent diagrams, together with the 
Gaussian decay ansatz (solid line Qcauss+NL)- We also plot 
for comparison the standard one-loop results, from the ratio 
(|139p (lower dashed line Qi), and from its expansion (|140p 
(upper dashed line Q2), as well as the phenomenological 
model of Pan et al. (2007) (dotted line Qpcs)- 



T^^'j'^''^'^\LUjD). In this fashion, we only need the nonlin- 
ear matter power-spectrum given by numerical simulations 
at the redshift of interest, where we wish to compute the 
equal-time bispectrum. Then, the equal-time reduced bis- 
pectrum Q is computed by using also the fit Pnl of the 
nonlinear power spectrum (from Smith et al. 2003) in the 
denominator of Eq. (|137p . We could apply the same proce- 
dure using the cosine time-dependent cutoff of Eas. (|109p . 
(|117[) , however since it does not fare as well as the Gaussian 
cutoff, we only show in Figs. [^01 1 ^ our results obtained with 
the Gaussian decay. 

In principle, we could improve this procedure in two 
ways. First, we could use in the external lines the actual 
two-time density correlation measured in simulations, in- 
stead of the Gaussian ansatz. Second, we could also use 
this in the internal lines. However, we leave this for future 
works, since the different-times density correlation has not 
been studied in details in current numerical simulations, 
and this would imply heavier computations for the relevant 
diagrams as explained above. Nevertheless, note that the 
different-time Gaussian decay agrees reasonably well with 
numerical simulations for the response function (Crocce & 
Scoccimarro 2007b) and it should also provide a reasonable 
guess for the nonlinear correlation function as discussed in 
Valageas (2007b), since it describes the random advection 
of small-scale structures by the larger-scale coherent flow 
which affects both two-point functions. 

We show our results for equilateral triangles as a func- 
tion of wavenumber k in Fig. [201 (curve Qcauss+NL)- For 
comparison, we also display the predictions of the standard 
one-loop perturbation theory (curves Qi a-nd Q2 associated 
with Eas. (|139p - p40p ) and the phenomenological model of 
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Fig. 21. The reduced bispectrum Q as a function of the 
angle 6 between wavenumbers ki and k2, as in Fig. [TOl 
but using the nonlinear power spectrum from Smith et al. 
(2003) for the two-point correlations within the one-loop 
steepest-descent diagrams, together with the Gaussian de- 
cay ansatz (sohd hue Qcauss+NL)- 



Pan et al. (2007). We can see that using the "exact" nonlin- 
ear power spectrum Pnl does not change much the results 
as compared with Fig. 1181 The failure to match the numer- 
ical simulations beyond 0.2h Mpc~^ shows that in order to 
improve this prediction it is not sufficient to improve the 
two-point functions alone. One needs to include additional 
diagrams to those displayed in Fig. (TUl such as those ob- 
tained by adding new lines that connect the external to the 
internal lines and which cannot be factorized in the form 
of Fig. [TO] (this requires new two-loop or higher-order dia- 
grams) . 

The dependence on the angle 9 between wavenumbers 
ki and k2 is displayed in Fig. [21] We only show the standard 
one-loop prediction Qi for comparison, because it fares bet- 
ter than the alternative Q2, see Fig. \T7\ Here we can see a 
significant improvement over the previous predictions since 
we partly recover the weakening of the dependence on 6 in 
the nonlinear regime. The improvement is most significant 
in panel (c) (where A|^(fc) = 1.3). At smaller scales the 
prediction Qcauss is too low as compared with the simula- 
tions, since it decays as in Eg. pSOp . as seen in Fig. [^Dj too. 
Thus, although it did not make much change for equilateral 
configurations fFig.l20p. using accurate two-point functions 
within such frameworks can improve the detailed angular 
behavior of the bispectrum in the weakly nonlinear regime. 

This means that resummation schemes can indeed im- 
prove over the standard perturbation theory (as shown by 
the comparison with the standard one-loop prediction Qi in 
Fig. [H}. However, phenomenological models (such as Pan 
et al. 2007) can still offer a better match to numerical sim- 
ulations. This is not completely surprising since the resum- 
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mation schemes discussed in this paper are still cxpectcid to 
break down at small scales, would it be only for the effect of 
shell-crossing that invalidates the equations of motion from 
which they arc derived (although one could apply similar 
methods to the Vlasov equation, sec Valageas 2004). On 
the other hand, being derived in a systematic fashion from 
the equations of motion, perturbative methods (that can 
involve various resummation procedures) are better con- 
trolled. Therefore, they should be useful and more reliable 
than simple models in their regime of validity. 

9. Conclusion 

In this article, we have described several methods to ob- 
tain systematic expansions for two-point and three-point 
correlation functions for the cosmological large-scale den- 
sity field. These methods also extend to higher-order cor- 
relations and response functions (e.g., the response of the 
two-point and three-point correlations to external pertur- 
bations or to the initial conditions). They also apply to 
related systems, defined by similar hydrodynamical equa- 
tions, such as the Zcldovich dynamics. We showed how the 
simplest perturbative expansion derived from the the path- 
integral formulation of the system (i.e. expanding over the 
cubic part of the action) corresponds to the equivalent stan- 
dard expansion directly obtained from the hydrodynami- 
cal equations of motion. Next, we described two rcsumma- 
tions, associated with "large-A^" limits, that can be conve- 
niently derived from the path integral (Valageas 2007a,b), 
and a high-fc resummation that is most easily derived from 
the equation of motion (Crocce & Scoccimarro 2006a, b; 
Valageas 2007b). Thus, we showed how to derive expansions 
for high-order correlations (or responses) and we compared 
the resummations involved in these various schemes. 

We first applied these methods up to one-loop order to 
the simpler case of the Zeldovich dynamics, where explicit 
expressions are simple enough to be given here. We explic- 
itly checked that the perturbative expansion of the path 
integral recovers the standard perturbative expansion over 
powers of the linear density field. However, since these two 
expansions are organized in different manners, we found 
that two diagrams, obtained within the former method for 
the three-point correlation, exhibit an additional UV di- 
vergence (for ri(oo) > —3). This divergence cancels out as 
we sum over all diagrams but it shows explicitly that some 
care must be taken in using these expansions. For numerical 
purposes, it is best to first regularize each diagram by sub- 
tracting the divergent parts that cancel out. Next, we ap- 
plied the large- A'' steepest-descent method to this Zeldovich 
dynamics. Wc showed how the time-dependence can be fac- 
torized, through functions T{u)jD), that measure the decay 
of different-time correlation and response functions in the 
nonlinear regime (within simple approximations). We found 
that, contrary to the standard expansions, tree diagrams 
and one-loop diagrams obey the same scaling laws over the 
linear amplitude Dag. This is due to the cutoff T{ujD) as- 
sociated with the partial resummations. We also explained 
how it is again possible to regularize the expansion, as for 
the previous expansion. 

Second, we applied these expansion schemes to the 
gravitational dynamics. All properties obtained for the 
Zeldovich dynamics remain valid and we have compared 
our results with numerical simulations. We checked that 
on quasi-linear scales all expansion schemes agree with one 



another and with the numerical simulations, up to 0.2h 
Mpc~^ (where A|(fc) ~ 1 for a ACDM power spectrum 
with fTg = 0.9). At smaller scales they deviate from one 
another and from the numerical simulations. Wc recovered 
the known failure of the standard perturbative expansion to 
reproduce the observed flattening in the nonlinear regime 
of the reduced bispectrum Q{ki, k2,k^) as a function of the 
angle between wavenumbers ki and k2. We found that the 
steepest-descent resummation does not fare much better. 
Even substituting the Gaussian decay obtained in a high- 
k limit for the different-time two-point functions does not 
help (it actually increases somewhat further the variation 
with 9). Nevertheless, the agreement with numerical simu- 
lations significantly improves if we use the "exact" nonlin- 
ear power spectrum (obtained from a fit to simulations) for 
the two-point correlation within this formalism (together 
with the Gaussian decay). This does not extend the range 
of validity for equilateral configurations towards higher k, 
but it manages to partly recover the observed flattening of 
the reduced bispectrum as a function of 9. 

We have also compared the predictions obtained from 
these expansion schemes with the simple phenomenological 
model of Pan et al. (2007). It appears that the validity of 
the latter extends to smaller scales, up to 0.8h Mpc~^ where 
A|(A;) ~ 17 in the case studied here. Since this model is also 
much simpler to use, it seems that the scale transformations 
devised by Hamilton et al. (1991) may indeed provide re- 
markably efficient tools to investigate various statistics as- 
sociated with gravitational clustering. At quasi-linear scales 
k ~ O.lh Mpc~^, the systematic expansion schemes stud- 
ied in this article may be somewhat more accurate (since 
they arise from rigorous systematic expansions) whereas the 
phenomenological model of Pan et al. (2007) may lead to a 
slight overestimate of the bispectrum (which would not be 
surprising since its expansion is not expected to agree with 
perturbation theory). 

The results obtained in this article suggest that, in order 
to make further progress along these lines, one needs to in- 
clude higher-order diagrams (two-loop diagrams or higher) 
to obtain reliable predictions for the bispectrum that ac- 
curately reproduce the weakening of the dependence on 
the angle 9 in the nonlinear regime. Moreover, using ac- 
curate two-point functions provides a significant improve- 
ment. This requires both the use of a Gaussian decay for 
different-time quantities (or a similar cutoff that agrees 
with numerical simulations) and of a reasonable equal-time 
correlation (e.g. using the nonlinear power spectrum mea- 
sured in simulations). On the other hand, the path-integral 
formulation might serve as a basis for other approximation 
methods than the expansion schemes discussed here. 

Alternative methods, based on phenomenological mod- 
els, such as the halo model (Cooray & Sheth 2002) or the 
scaling transformation of Hamilton et al. (1991), may pro- 
vide other means to describe the statistical properties of 
gravitational clustering. Although simpler and more effi- 
cient in the nonlinear regime, they do not have the rigorous 
and systematic character of expansion schemes. Moreover, 
they often involve several free parameters that must be fit- 
ted to simulations and they are not derived from first prin- 
ciples. It would be interesting to develop a connection be- 
tween such phenomenological approaches and more system- 
atic methods as those developed in this article. However, it 
is not clear whether such a connection can be achieved at 
a detailed level. This issue is left for future works. 
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Appendix A: Bispectrum for the gravitational 
dynamics: tree and one-loop diagrams for the 
expansion over Kg 

We give in this Appendix the expressions of the diagrams of 
Fig. [HI associated with the expansion over powers of Kg of 
the path integral defining the bispectrum, for the case of the 
gravitational dynamics. The corresponding results for the 
Zeldovich dynamics are given in Sect. [721 The expressions 
below are similar to the former but the rational functions 
which appear in the integrands are more intricate. First, 
the tree-diagram (a) is: 



B('^)(fcl,A;2,fc3) =-PL(fc2)-PL(fc3) 
, 4(k2.k3) 



"10 


(k2 




k2.k3 


-+ 




k2j 


k2k3 



Ik^k"^ 



2 perm. 



The diagram (f) which involves a triple integration over the 
linear power spectrum reads as 

B'^^'^ = y"dqidq2dq3 ^^(qi + q2 - ki) 





f q3 + 


k2)PL{qi)PL{q2)PL{q3) 


IQqlql 




+ 9i)(qi-q2) - 


M(qi.q2)2 












^nqi 


+ g|)(qi-q3) - 


M(qi.q3)2 












-nqi 


+ 9i)(q2-q3) - 


|-4(q2.q3)2 





(A.2) 



Next, the two diagrams (b) and (g) which involve a double 
integration over are: 

S^*") - -PL(fcl) y"dqidq2 <5D(qi +q2 -k2) 
lOqfqj + Tjqf + gi)(qi.q2) + 4(qi.q2)' 

7A:|(ki.k3)[10<725i + 7(g2 + '?i)(qi-q2) + 4(qi.q2)2] 



+ 



25 2 



2b2klklqlql 
6- 



7fc2(k2.k3)-2fc|(ki.k2) L , 7{ql + ql){ci,.Ci2) 



9 9 
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9 9 
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PL{qi)PL{q2) + Sperm. 



(A.3) 



B^^^ = Pdki) J dqidq2 Soim + q2 - k2) 

10(k2.qi)(k2.q2) -3A:2(qi.q2) 
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9 9 
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882fc2g2 
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q2| 
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-Sperm. (A. 4) 



Finally, the five diagrams (c), (d), (e), (h), and (i), which 
involve a single integration over P^, are given by: 



R(c) p 10fc|fc2-f7(fc|-ffcg)(k2.k3)+4(k2.k3)2 
'=PL[k2)PL[k3) — jp 

' '^2 "-3 

"30A:2g2 + (47fc2_7g2)(j^2.q) 



X dqPL{q) 211 12 

J [ 126(72|k2 - q|2 
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(A.5) 



(A.l) B^'^' =PL(fc2)PL(fc3) / dqPi(g) 
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"2(k2.q)^(4(k2.q)-5g^)-ffc|(6g^-29g^(k2.q) + 21(k2.q)^) Because of the property ki + ka + kg = 0, which has been 

used to simplify somewhat these expressions, it is possible 
to write each contribution under many different forms (for 
instance by replacing ki by — (ka -f kg)). 

We do not give in this paper the expressions obtained 
within the steepest-descent expansion, as they are longer 
than Eqs. (jA.l|) - (jA.9p . However, as explained in the main 

text, by putting all time-dependent functions T^^^^^ equal to 
unity, they reduce to the expressions above for diagrams 
(a), (f), (g), (h), and (i). 
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g2|k2-q|2 




6fc3V + 7(A:|+q2)(k3.q) 


+ 8(k3.q)2\ 



92|k3-Hqp 

B^'^=PL{k2)PL{k3)fdqPL{q) 



ki.q 



-5 perm., (A. 8) 



2fc|g^+(fci + g^)(k2.q) 

2„2 



kjq 



2k3 . (k2 + q) |ki - q| ^ - 7kl [ki .k2 + q. (ki - k2 - q)] 
154fc||k2 + qp 
5k3.(ki-q)\ , 2k2.q / k3.(ki ~q)|k2+qp 



66k^q'^ 



2k3.(k2+q)|ki-q|2 - 7fc|[ki.k2 + q.(ki-k2-q)] 



231fc|Q2 



+ 



7fc,V + (2fci^-9g^)(ki.q) 
g2|ki-q|2 



2fc2V + (fci + g^)(k2.q) 

fc2g2 



•2 k3.(k2 + q)|ki-qp - fc|[ki.k2 + q.(ki-k2-q)] 
462fc||k2 + q|2 
5k3.(ki-q)\ , 2k2.q / k3.(ki -q)|k2+q|2 



1386fc^ 



^2 



1386/c|g2 



2k3.(k2 + q)|ki-q|2 - /c|[ki.k2 + q.(ki-k2-q) 



693A:|g2 



-1-5 perm. 



(A.9) 



